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Abstract. This is the first of a series of three papers treating light baryon resonances (up to 3 GeV) 
' within a relativistically covariant quark model based on the three-fermion Bethe-Salpeter equation with 

^Sj , instantaneous two- and three-body forces. In this paper we give a unified description of the theoretical 

background and demonstrate how to solve the Bethe-Salpeter equation by a reduction to the Salpeter 
equation. The specific new features of our covariant Salpeter model with respect to the usual nonrelativistic 
^ , quark model are discussed in detail. The purely theoretical results obtained in this naper will be applied 

■ numerically to explicit quark models for light baryons in two subsequent papers jl|,|^. 
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, 1 Introduction 

1^ _ The classification of baryon resonances as three-quark states within nonrelativistic potential models has a long and very 
successful history. It is however unclear how to relate such models to QCD. Some ingredients of nonrelativistic quark 
models emerge from QCD, e.g. massive quarks as a consequence of chiral symmetry breaking, linear confinement 
potentials (on the lattice) due to the nonabelian gauge coupling and some candidates for spin-dependent residual 
interactions like one-gluon-exchange or instanton-induced quark forces. For light quark flavors it is however unclear, 
k>( ■ how to unite these features in a common picture. The main obstacle is the nonrelativistic approach which seems to be 
■ completely inadequate for small constituent quark masses and strong quark binding. 

Quantum field theory seems to offer a solution to this problem, replacing the nonrelativistic wavcfunctions by 
Bethc-Salpeter amplitudes obeying a suitable Bethe-Salpeter equation. In the case of QCD none of the basic ingredi- 
ents of these equations is reliably known, i.e. we have no reliable prescription to calculate the full quark propagators 
and interaction vertices. Moreover we meet a serious problem with gauge invariance because the Bethe-Salpeter am- 
plitudes are gauge-dependent. Nonetheless the general framework of quantum field theory can be used for a reasonable 
phenomenological description. If we want to remain as close as possible to the features of nonrelativistic quark models 
the Bethe-Salpeter equation should contain free quark propagators with constituent quark masses and instantaneous, 
unretarded interactions only. Both requirements are purely phenomenological assumptions but reasonably justified by 
the apparent success of nonrelativistic quark models. In this way these Bethe-Salpeter amplitudes form a more suit- 
able basis for quark models, but respecting, in particular, relativistic covariance. As such it was already successfully 
used for the description of light mesons |jl^,|l9,^|2^,^,^|2^]. The baryon Bethe-Salpeter equation with genuine 
instantaneous three-quark forces is solved as in the mesonic calculations by a reduction to a three-dimensional integral 
equation (Salpeter equation) which is very similar to the Schrodinger equation. The spectrum contains however also 
antiparticle solutions corresponding to particles with charge conjugated quantum numbers. This situation is new and 
needs a special discussion. Another complication arises when genuine two-particle interactions are taken into account. 
In quark models this is natural, when the (three-body) confinement forces are supplemented by a two-body residual 
interaction (one-gluon-exchange, instanton induced forces). In this case an effective three-body interaction kernel has 
to be derived. 

None of these features is entirely new, but there is no reference in the literature which presents this theoretical 
background in a unified way. The purpose of this paper is to fill this gap. In two consecutive papers we will use 
these purely theoretical results for specific calculations of the baryon spectrum up to 3 GeV. 
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This paper is organized as follows: In section]^ we briefly recall how in quantum field theory bound states of three 
fermions occur as poles in the six-point Green's function defining the Bethe-Salpeter amplitudes as the corresponding 
residua at these poles. This property of the Green's function is used in section ^ do derive simultaneously the Bethe- 
Salpeter equation for the Bethe-Salpeter amplitudes and their normalization condition in a simple and appealing way 
by a Laurent expansion of the integral equation for the six-point Green's function in the vicinity of this pole. Section 
^ is concerned with the reduction of the full eight-dimensional Bethe-Salpeter equation to a six-dimensional Salpeter 
equation by integrating out the relative energy dependence of the full Bethe-Salpeter amplitudes. To this end we use 
a covariant formulation of the instantaneous approximation for three- and two-body interaction kernels and assume 
that the full quark propagators can be suitably approximated by their free forms introducing effective constituent 
quark masses. In a first step, taking only the genuine (instantaneous) three-body kernels into account, we show how 
a straightforward reduction can then be performed, thus yielding a reduced equation which may be formulated as 
an ordinary eigenvalue problem in Hamiltonian form, where the Hamiltonian is hermitean with respect to a scalar 
product induced by the normalization condition of the Salpeter amplitudes. Complications arise for the more general 
case when also genuine two-particle interactions are taken into account. This case needs a special discussion and we 
demonstrate that a reduction to a Salpeter equation in the same Hamiltonian form can nevertheless be achieved by 
deriving an effective instantaneous three-body kernel which parameterizes all retardation effects of the unconnected 
two-body interactions. In section || we present the Salpeter equation in Born approximation of the quasi potential 
which constitutes the basic covariant equation of our model. We discuss the structure and main features of the Salpeter 
equation and its solutions with respect to the ordinary nonrclativistic quark model. Special features discussed in this 
section are the one-to-one correspondence of the Salpeter amplitudes to the states of the nonrclativistic quark model 
and the additional anti-particle solutions of the Salpeter equation. Finally we give a summary and conclusion in section 



2 Green's functions and Bethe-Salpeter amplitudes 

In non-relativistic quantum mechanics, a bound state of three particles is described by a normalized wave function 
satisfying the three-body Schrodinger equation. This is in general the underlying equation for the description of baryons 
as bound states of three quarks in the framework of the various phenomenological non-relativistic potential models. A 
more profound basis for describing bound states in relativistic quantum field theory is the Bethe-Salpeter equation |^ 
for the so-called Bethc-Salpeter amplitudes, which might be considered as the covariant analogues of 'wave functions' 
in the non-relativistic case. The Bethe Salpeter equation has been first derived for the two-particle system by Salpeter 
and Bethe Taylor Q investigated the application of the Bethe-Salpeter equation to the three-body system. 

In this section, we outline a method to treat the three fermion bound state problem in relativistic quantum field 
theory by using Green's function techniques. This allows to derive the Bethe-Salpeter equation for three bound fermions 
simultaneously with the normalization condition of the corresponding Bethe-Salpeter amplitudes. The method is based 
on the fact that in general a bound state of elementary particles, whose fields appear in the underlying interaction 
Lagrangian, corresponds to a pole in the total energy of the Feynman propagator (Green's function) of the many 
particle system. These poles do not arise from single perturbative Feynman diagrams, but rather from an infinite 
series of diagrams. In this context the Bethe-Salpeter amplitude is then defined as the residuum of the bound-state 
pole of the Green's function. This connection between bound states and the singularities of Green's functions was 
originally the basis of the first rigorous proof of the two-particle Bethe-Salpeter equation given by Gell-Mann and Low 
1^. However, this non-perturbative approach is clearly general and can be applied generically to the n-body Green's 
function as shown for instance in the textbook of Weinberg As mentioned above, we apply this method to the case 
of three fermions (quarks) only. It consists of the following three steps ||^,^,^: 

1. Sta rting point is the six-point Green's function describing the propagation of three interacting fermions. In section 



2.1 



we analyze the structure of the usual perturbative power series expansion of the three-quark Feynman propaga- 
tor: introducing the concept of irreducible interaction kernels for the case of three particles in a manner similar to 
that of Salpeter and Bethe in the two-particle case ||^, we outline, how the infinite power series can be rearranged 
into an inhomogeneous integral equation. 



In section 2.5 we examine the analytical structure of the six-point Green's function: we isolate the contribution 
of a three-fermion bound state to the six-point Green's function and show, how the bound state gives rise to a 
pole in the total energy variable (or in the invariant total four-momentum squared). This procedure defines the 
Bethe-Salpeter amplitudes of a specific bound state by the residue of the corresponding bound-state pole which 
factorizes into the Bethe-Salpeter amplitude and its adjoint. 

Finally, by a La uren t ex pansion of the Green's function in the vicinity of this bound-state pole and using the results 



from sects. 2.1 and 2.5, we will derive a homogeneous integral equation for bound states, i.e. the Bethe-Salpeter 
equation along with the normalization condition of the corresponding amplitudes. This will be done in section ^. 
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2.1 The six-point Green's function for three fermions 

The fundamental quantity describing three interacting fermions in quantum field theory is the six-point Greens's 
function (or three-fermion Feynman propagator), which is the vacuum expectation value of a time ordered product of 

three fermion field operators and their adjoints ^ := S/'^'yO in the Heisenberg picture: 

Craia2a3; a[a'^a'^{xi, X2, X3; Xi, X2, X^) := (1) 

~(0| T <(xi)^'f,(a;2)ff3(x3)^^^K)^^, (x^)^:;, (4)|0). 

Here the ~ {ai,fi,Ci) denote multi-indices combining the indices of the quark fields ai in Dirac, fi in flavor and 
Ci in color space. |0) denotes the true physical vacuum state and T is the time ordering operator acting on a general 
n-fold product of Heisenberg fermion field operators oi 'F, {i = 1, . . . ,n) defined as 

T{A\xi)A^{x2)---A"{x^)} 
= sign(a) T { ) A'^^^) ^^^^^^ ) . . . ^.(n) ^^^^^^ ^ | 

where the sum runs over all permutations a G Sn with signum sign(cr). is a generalization of the usual Heaviside 
function 

n n n ( 1 iOT X^ > X^ > ■ ■ ■ > X^ 

0{x'„xl...,x^,) = \ ^ • (3) 

I otherwise 

Among other possibilities (depending on the time ordering considered) the six-point Green's function G represents 
the probability amplitude for three (generally off-shell) quarks to propagate from space-time points x'l ,X2,x'^ to 
Xi^X2,X3. Using the technique of ordinary time-dependent perturbation theory, the six-point Green's function G may 
be expressed in the form of an infinite power series (see any standard textbook of quantum field theory, for instance 
0): 



G{xi,X2,X3;x[,X2,x'3) ^ — : r y2—n- / d'*?/! . . . d^2/fc 

"^\Te^p(~if^^dtHj,it))\0)t[ ^' 



x{0\TW}^{x,)^lix2)^f^ix3)^]pix[)^^^^^ . ..nM\0). (4) 

Now the state |0) represents the unperturbed vacuum and !f/p, tf^/p, Hjp and Tiip are the field operators, the interaction 
Hamiltonian and the Hamiltonian density operator in the interaction picture, respectively. 

2.2 The integral equation for the six-point Green's function 

Using Wick's theorem for time ordered products of field operators, the right hand side of eq. (^) may be evaluated 
order by order (in the coupling constant) to obtain a power series expansion which may be represented in terms of 
ordinary Feynman graphs describing the interaction of two or three fermions in finite order (see fig. |l]) . 




+ + + 



Fig. 1. Finite order perturbative contributions to the six-point Green's function G. 
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In scattering processes (at high energy), where neither a three-body bound state nor a two-body bound state 
in any of the two-particle subsystems occurs, only a finite set of diagrams may be taken into consideration. The 
investigation of bound states, however, requires to go beyond such a perturbative approach, i.e. an infinite sum of 
diagrams (or at least an infinite subset of diagrams) has to be taken into account. The reason for this is that e.g. a 
thr ee-body bound state leads to a pole of the Green's function in the total energy variable, as we will see in sect. 
2.5. But such a pole never arises from a finite set of Feynman diagrams alone. To go beyond perturbation theory, 
one recasts the infinitepower series expansion (^ in the form of an inhomogeneous integral equation, as it was done 
by Bethe and Salpeter 0] for the case of two particles. Let us briefiy sketch this procedure for the case of three fermions: 



1. ) One introduces the concept of irreducibility, i.e. one classifies all those diagrams appearing in the power 
expansion series (^ in reducible and irreducible graphs. For the definition of (ir)reducibility in the case of three 
interacting particles we distinguish two- and three-particle interactions: 

— A connected two-fcrmion interaction graph is called irreducible, if it cannot be split into two sirnpler graphs by 
cutting two fermion lines only. Some examples of irreducible two-body diagrams are shown in fig. ^. 

— Correspondingly, a connected three-fermion interaction is called irreducible, if it cannot be separated into two 
simpler graphs by just cutting three fermion lines. Examples of such graphs are given in fig. ^ 

— All other interaction graphs are called reducible. Clearly, due to the above definitions of irreducibility, reducible 
diagrams can always be cut into irreducible parts. 

2. ) The (infinite) sum of all irreducible connected two-particle graphs is collected into the so-called irreducible 
two-particle interaction kernel 

(5) 



K 



aia2; a^aJj 

See fig. ^ for a diagrammatic representation of K^'^"' . Similarly, all irreducible connected three-particle graphs are added 
up to the so-called irreducible three-particle interaction kernel 



K 



A graphical picture of K^^"' is shown in fig. ^. The arguments x'^, xi and multi- indices 
the coordinate space and the Dirac, flavor and color space dependences of the kernels, respectively. 



a,-, ai m 



(6) 

@) and (H) indicate 




+ 



+ 



+ 



+ 



Fig. 2. Graphical representation of the two-particle irreducible Bethe-Salpeter kernel K'^'^^ as sum of all possible connected 
irreducible two-particle interactions. 




+ 



+ 



+ 



+ 



Fig. 3. Diagrammatic picture of the three-particle irreducible Bethe-Salpeter kernel iC'^-* as sum of all possible connected 
irreducible three-particle interactions. 



3.) Apart from the connected two- and three-particle interactions, applying Wick's theorem to the right hand side 
of eq. (^) also generates unconnected terms, as e.g. the bare quark propagators, but moreover all kinds of self-energy 
contributions to the single fermion lines of each quark, summing up to the full quark propagators 

S'Fa.a'S^^^<) = {0\T^:^ix,)K'ix'm , (7) 

as indicated in fig. 0. 
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Fig. 4. Perturbation series of the full dressed quark propagators defined in eq. 



4.) All reducible interaction diagrams of any desired order appearing in the power series expansion can now be 
generated by iteration of the irreducible two-particle (in each quark pair) and three-particle interaction kernels if^^-* 
and using the full quark propagators Sp for the inner fermion lines. This is accomplished to all orders by virtue 
of the following inhomogeneous integral equation, which uses the two- and three-particle interaction kernels as integral 
kernels [|,|], i.e. 

^010203; a'^ai^a'^ •^2i ^3; -^l? -^2 7 -^s) 



i J d'^yi d^y2 d^ys S"];. aibi (a^i, yi) S'l 0262 (^^2, ^2) S], ^^f,^{x3,y3) 
J dVi d^'ya d^'yg Klf^^^^.^,^^,^^,^{yi,y2,y3;y'i,y'2,y3) Gv^v^u^, a',a'^a'^{y\,y'2.y'?,-,x\,x'2,x'3) 

^ / dV d''y2 S], ^^f,^{xi,yi) S% ^^t,^{x2,y2) 

cycl.P„rm. 

4„/ ^4 / 7^(2) / I i\ ^ 11 I ^ . ^1 I I- 



-1 

cy 

(123) 



d 2/1 d 2/2 -^^^6162; b;^(yi'2^2;yi,y2) Gb;^a3; a;4a!,(yi,2/2:a;3;2:'i,4>4)> (8) 



where our notation implies summation over indices bi, b[ occurring twice. For a diagrammatic illustration of this 
integral equation, see fig. |[ In fact, the Neumann iteration of this integral equation reproduces all possible reducible 
interactions and thus precisely all the terms of the power series expansion of eq. (B). 



G 




cycl. perm. 

X,, X,, X, 









X, 






Xi 








< 




Fig. 5. Graphical illustration of the inhomogeneous integral equation ^ for the six-point Green's function G. JC'^-* and K''^^ 
denote the irreducible three- and two-body Bethe-Salpeter kernels, respectively, represented graphically in figs. ^ and ^. Thick 
arrows on quark lines indicate full dressed quark propagators as shown diagrammatically in fig. ^ 



Note that also the irreducible interaction kernels K'^'^^ and K'^^'^ consist already of an infinite number of graphs 
and in general cannot be calculated exactly. They are basically unknown functions and thus have to be parameterized 
phenomenologically. However, the decisive advantage of the non-perturbative construction of the Green's function G 
from an inhomogeneous integral equation (|^) is that its solution automatically implies an infinite number of interac- 
tions even if the kernels are approximated by their lowest order Born terms, which constitutes the so-called ladder 
approximation. Such an approximation is sufficient in theories, where the coupling constant is small and the interaction 
kernels may be considered as an asymptotic series expanded in terms of the (small) coupling constant. (In such a case 
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one would expect most of the binding of a bound state to come from the repeated action of the Born diagrams alone.) 

For further discussion of eq. it is useful to introduce an appropriate compact notation. First let us combine the 
irreducible two- and three-body kernels if and K^'^'' to a single integral kernel K. We introduce the inverse Sp ^ 
of the full quark propagator Sp by 

j dV ^Fa,b(a;fc,yfc) S^;^,^{yk,x'k)^Sa,a'^ (5(4)(Xfc-4). (9) 

This allows to rewrite the sum of the two-particle interactions ii'^^^ in each quark pair in the form of a three-body 
kernel 

(ijk) = 

cycl. perm, of 
(123) 

In this form we can combine the two-body interaction kernels with the three-body kernel K^^"^ to a uniform integral 
kernel K (see fig. ||): 



iK = 



iK 



+ E 



(2)- 



cycl. perm. 



=: iK' 



(2) 



Fig. 6. The integral kernel K combining the three-body irreducible kernel and the two-body irreducible kernels in each quark 
pair; the filled circle denotes an inverse full quark propagator. 



Moreover, we introduce the symbol Go for the triple tensor product of the single quark propagators Sp which is the 
lowest order contribution to G: 

Go a,a2a,:a',a'y^ixi,X2,X3;x[,x'2,x'3) := S], ^^^,^{xi, x'^) Sj, ^^^,^{X2, x'^) S^^ ^^^,^ix3 , x'^) . (12) 

Finally, we define a shorthand operator product notation for the summation over indices and the integral operation 
in coordinate space: 

-^]aia2a3; a'^ajjof, ^aia2a3; 616263 -^616263; o'lO^a!, i (13) 

[A B]{xi,X2,X3;x[,x'2,X3) J d^yi(fy2d'^y3 A{xi,X2,X3;yi,y2,y3) B{yi,y2,y3;x[,x'2,X3). (14) 

With these definitions the inhomogeneous integral equation for the six-point Green's function can be represented in 
the more compact form of an operator equation 

G = Go-iGoKG. (15) 

Note that this integral equation for the Green's function G can also be written in its equivalent adjoint form, where 
the operator product Go if G on the right hand side of eq. ( p^ ) appears in reverse order: 

G-Go-iGifGo. (16) 

The equivalence of the integral equation ( |l5| ) and its adjoint ( |l^ ) is obvious, since both equations have the same 
Neumann series. 
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2.3 Space-time translational invariance 



The six-point Green's function G as defined in eq. is invariant under arbitrary space-time translations, i.e. 
G{xi , X2, X3; x[, X2, x'^) = G{xi + a,X2 + a, X3 + a; x[ + a,X2 + a, x'^ + a) 



(17) 



for all a £ 11"*. Due to this symmetry it is natural to introduce new coordinates, namely an external 'center-of-mass' 
coordinate X and internal, i.e. translationally invariant, relative coordinates ^ and ry, the so-called Jacobi coordinates. 
We choose: 

X := UX1+X2+X3), Xi=X + ^^-'^" 



Xl - X2, 



7] := +X2- 2x3), 



X2 



X3 



(18) 



X 



The space-like components X, ^ and rj of these variables can be interpreted in the non-relativistic limit as usual center- 
of-mass and relative coordinates for a system of three particles with equal mass. However, in a covariant framework 
this choice is a priori arbitrary and the variables X, ^ and rj have in general no direct physical meaning. Choosing 
now specifically a :— —^{X + X') in eq. (|l^ ) we find that in fact the six-point Green's function G depends only on 
translationally invariant coordinate differences X — X', ^, r/, ^' and 77', i.e. 



G{xi,X2,X3; x[,X2,X3) = G{X - X';^,ri;C,v')- 



(19) 



Of course, the same holds also for the triple product Gq of the free single quark propagators, and the translation 

invariance of the Green's function G necessarily implies that in particular the interaction kernels K^^") and K 
must by themselves be translationally invariant quantities. In momentum space space-time translation invariance is 
equivalent to the conservation of the total four-momentum. Consequently, as will be shown in the following subsection, 
the twelve-dimensional integral equations ( 15) for the six-point Green's function in coordinate space and its adjoint ( |l^ ) 
after Fourier transformation become only eight-dimensional integral equations in the momentum space representation. 
Due to momentum conservation, these momentum space representations depend only parametrically on the total four- 
momentum. To perform the Fourier transformation let us define the corresponding conjugate momenta to X, ^ and rj 
which are given by the total four-momentum P and the following relative four- momenta and Pr/'. 



The new sets of coordinates 



P pi +P2 +P3, 

Pi ■= \{pi -P2), 

Pv •= ibi +P2- 2^3), 



Pi ^ ^P + Pi^ 

^ P2 = \P-Pi^ 
P3 = ^P-p,,. 



and (EG) satisfy the condition 



(20) 



{pi,Xi) + {P2,X2) + {P3,X3) {P,X) + (p^,0 + {Pn,V), 



(21) 



and a technical advantage of this special choice of variables is that the Jacobians of the transformations (|l^) and (|2 
are unity, i.e. 



d{X,^,v) 



d{xi,X2,X3) 



1 and 



diP^P^^P-n) 



dipi,P2,P3) 



1. 



(22) 



2.4 Momentum space representation of the integral equation 

For any six-point function A = G, Gq, K, K and K^'^\ i.e. the six-point Green's function, the triple product of 
quark propagators or the interaction kernels, we define the Fourier transform by 

[TA] {pi,P2,P3;p'l,P2,P3) ■ = 

d^xi d*X2 d^X3 e+i«Pi^^i>+<P='^^>+<P^'^^» (23) 
d^x[ d^x'2 d^x'3 e-H<Pi.<)+(pi.4>+(P3,4>) A{xi,X2,X3;x[,x'2,x'3). 
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Using the properties ( pT| ) and ( p2| ) of the new coordinate sets, the Fourier transforms can be written in terms of relative 
Jacobi momenta and the total four-momenta, 

[TA]{pi,p2,P3;p[,p'^,p'3) ^ Apip^,p^;p'^,p'^) (27r)V4)(P - P')- (24) 

Due to the translational invariance of the six-point functions A the (5-function reflects the conservation P' — P oi 
total four- momentum. The remaining part Ap, just depending parametrically on P, is defined by the following Fourier 
transformation 

A{xi,X2,X3;x\,x'2,x'^) = A{X - X';^,r];C 

d^-P -i{P,x-x') f p-i(p,,«) .-iK,^) f J{p'a')JipW) 



_ 'J^J^. - J ^ ^ ' Apte,p,;pJ,K,). (25) 

which exhibits the exclusive dependence on the relative coordinates X — X' , ^, ry, ^' and 77'. The momentum space 
representation Gop of the quark propagators Go then reads explicitly 

Gop iPi,Pr,;p'^,p'^) = 5^ i^P+Pi + hPn) ® Si {^P-p^ + ^p,,) ® Si {\P-Pn) 

X (27r)4j(4)(p^_p,) (27r)4 5(4)(p,-p;), (26) 

where (due to translational invariance) the Fourier transforms of the single quark propagators are defined by 

SUx^.x\) = SU^, - x[) / ^ e"i<^----^> S^Pr)- (27) 

(2) 

For the sake of completeness we should also specify the explicit form of the Fourier transform K p of the two-particle 

term if defined in eq. (|l^). To this end we first have to define the Fourier transform of the two-particle interaction 
kernel K^'^'^ {xi, xj; x[,x'j). According to translational invariance, it is useful to introduce two-particle 'center of mass' 
and relative coordinates and S^k for each possible quark pair (ij), i.e. 

Xk hixi + Xj) 

^ for [ijk) = cycl. perm, of (123), (28) 

C/c • — Xi Xj 

as well as their corresponding conjugate variables, the total two-particle momenta Pfc and the relative momenta p^^. , 
i.e. 

Pk ^^Pr+Pj = ^P+Pr„ 

Prik — Up^+P] - '2Pk) (29) 

Pik '■— \iPi ~ Pi) for {ijk) = cycl. perm, of (123). 

Note that we have expressed Pk by the total three-particle momentum P and an additional variable Pri^ in order to 
relate the sets (P^ , ) of two-particle momenta to the set (^0|) of relative three-particle momenta (p^ , P?) ) = {p^^ , Pr^ ) in 
the case [ijk] = (123) and to the equivalent cyclically permuted sets {p^i^Pth) and {p^2yPv2) in the cases {ijk) = (231) 
and (312), respectively. The cyclically permuted sets of the relative momenta are obtained by linear transformations 
of the existing set (EG) according to 

When K'^'^^ depends on translationally invariant two-particle variables only, the Fourier transform of if^^^ is given as 
K^^\x,,xf,x[,xr)^K(^\Xk-X'^;^k,ek) (31) 

(27r)4' J {2n)^ ]{2n)^ ^P^P^.^PiJ- 

Using the definition ( pO[ ) of if^^'' and the definitions ( p5[ ) and ( pl| ) of the Fourier transforms of k''^^ and if^^^, we find 
:it 

aia2a3;a[a'^a'^{PiiPv'^ P'i^P'r) = (32) 
E ^;|U,.Ja.a,;.'.'(P«.'^'L) ^^a^^ (l^-^.J (2-)^ 5i'\p..~p'J. 

(123), (231), (312) 



-(2) 

the following explicit form for K p 
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where Sp {pk) is the momentum space representation of the inverse of the fuU quark propagator defined in eq. 
which obeys 

E ^Fa^b, (P>^) ^FbU (P'^) = ■ (33) 

bk 



With definition (g5|) of the Fourier transforms of G, Go and K, the properties ( plD and (22) of the Jacobi coordinates 
and the exphcit form ( ^ ) of Gqp, we are now in the position to write the inhomogeneous integral equation (|l5|) for 
the six-point Green's function G in its momentum space representation, 

GpiPi,Pr,\p'^,p',,) = (iP + + ^Pr,) ® Si (iP -p^ + ip^) ® S% (iP - pr,) 

X (27r)4 5W(p^_p,) {27r)U^'\p,,-p'^) 



+ S^p (iP + + hPv) ^1 (i^ - P« + iPv) ® -^l, (iP - p^) 

("^^ / 1^ ^ Xp(P5,P,;p^',p;') Gp(p^',p:;;p^,p;), (34) 



where we suppressed the dependences on the indices using the shorthand tensor notation and the definition (13) of 
the operator product. Note that, due to the conservation of the total four momentum, the inhomogeneous integral 
equation depends only parametrically on the total four momentum P, while the integral operation involves only the 
relative momenta pe and p^. Let us therefore introduce the momentum space representation of the operator product 
corresponding to dl4) as 

[Ap Bp] ip^,p,,;p'^,p'^) := j ^ ^ v4p(p^,p,;p^',p;') Pp(p^',p;';p^,p;) (35) 



which again allows to write the momentum space representations of the integral equation (|15|) and its adjoint ( [16D in 
a concise operator notation: 

Gp ^ Gqp - i Gqp Kp Gp, (36) 
Gp = Gop — i Gp Kp Gqp, (37) 

with the subscript P indicating the parametrical dependence on the total four momentum, which becomes important 
in the next two sections for the investigation of bound-state contributions to G. Accordingly, we will evaluate eqs. ( |3^ ) 
and (|37| ) at the positions P — P, where bound states with mass Af^ = occur, allowing the derivation of the bound 
state Bethe-Salpeter equation and the normalization condition of the corresponding amplitudes. But to this end we 
first have to know how the six-point Green's function behaves at these bound-state pole positions. This is the topic of 
the next subsection. 



2.5 Bound-state contributions - Bethe-Salpeter amplitudes 

In quantum field theory bound states are related to the occurrence of poles of Green's functions in the total energy 
variable P° or, equivalently, in the invariant four- momentum squared P^. Here we verify this statement for the case 
of the six-point Green's function Gp. 

In the following we consider bound states of three quarks with (positive) mass M and positive energy wp := 
VP^ + M^. The corresponding Fock states with total four momentum P = (wp, P) and mass P^ = Af^ are denoted 
by |P). They are eigenstates of the total four-momentum operator P — pi + P2 + Ps, i-s-. 

P\P)=P |P), (38) 

and are normalized covariantly according to 

(P|P') = (2^)^2tjp5(3)(p_p'). (39) 

The six-point Green's function (|^) in general describes all possible kinds of processes with three incoming and three 
outgoing fermions. The 'initial' and 'final' fermion or anti-fermion lines, however, are not yet fixed, until a particular 
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time-ordering has been chosen. Here we are interested in the extraction of 'baryon' contributions to G, i.e. real bound 
states of three quarks with positive energy that propagate forward in time. Therefore we shall consider those specific 
contributions to the six-point Green's function G which have the particular time orderings x1,X2,x'^ > xf , xf , x';^ , 
i.e. which contain 

e (mm{x'i,x°2,x'i) ~ max{x'°,x'tx'l)) = (l for 4, 4, > x'^,, x'l x'l (40) 
V ^ 1 otherwise. 

Isolating this part of the Green's function defined in eq. (|l|), we have 

Gaia2a3; a'^a'^a'^i^li X^] X^, X2, X^) = 



T 



{^l{x,)^lix2)Wlix,)} T {^\{x',)^\{x'2)€>y,)}\0) 
\^mm[Xi,X2,X3) - max(a: i,x 2,x ^jj 



other terms arising from different time-orderings. (41) 



Now we can evaluate that contribution to the Green's function which arises from three-quark bound states ( p8D with 
mass M, by inserting the complete set of the intermediate states |P) in between the two time-ordered products in the 
matrix element (|4l|): 

" / (2^f^ {o\T^l^ix,)^lix2)^lixs)\P) (P| r<,K)<,(4)^a^(4)|o) 

X I mm(x{^, X2, X3) ~ max(a; i,x 2,x 3) J 

-I- other terms. (42) 

Here 'other terms' now denotes the terms not only arising from other time-orderings, but also from other intermediate 
states. 

We define the Bethe-Salpeter amplitude Xp for the bound state |P) and its adjoint xp by the following 
transition amplitudes between the state \P) and the vacuum |0), 

XPa,a,aM,X2,Xs) {0\ T ^I^ixi)^l{x2)^l{x3) \P), (43) 

Xpa',ai<(4,4,4) (^1 T¥l,(x[WUx'2WUx',) |0), (44) 



which appear in the bound-state contribution (|4^) to the Green's function G. Due to translational invariance we can 
factorize the total momentum dependence of the Bethe-Salpeter amplitude Xp ^^id its adjoint Xp which contributes 
just by a trivial phase factor: 



Xp 



p{xi,X2,X3) = e '^-f''-^^ Xp{£„v) 



Xp[Xi,X2,x^) ^ e^ ' Xp[^,Vj 

(27r)4 (27r)4 

Thus, we obtain translationally invariant Bethe-Salpeter amplitudes and their Fourier transforms which depend only 
on the internal relative coordinates ^, r] and p^, Pr,, respectively. 



i(^.^'> f ^ ^ e'ip'vO e'H-V) xpivlp',). (46) 
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The 6*- function in eq. (^2|), which dictates the specific time ordering x1,X2tXt^ > x'i,X2,x'^, gives rise to a pole of 
G in the total energy variable P° and we finally arrive at the following Laurent expansion of Gp in momentum space 
near the pole at P° = uip: 



or written covariantly 



2wt 



Xp{p^,P-n) Xp{p'^,p'rj) 
po - wp + ie 



regular terms for — > cjp, 



Gpip,,p,,p',p',) ^ -i + regular terms for P^ 



(47) 



(48) 



where we have introduced a six-point function [xp XpI by the separable product of the Bethe-Salpeter amplitudes 
allowing us to suppress the dependence on indices in ( p7| ) and (^): 



XpiPi,Pr,) Xp(j>'(,,Pn) 



Xp aia2CL3 {Pi,Pv) Xp 



(49) 



This typical analytical structure of the six-point Green's function Gp in the vicinity of the bound-state P P 
enables us to isolate the three-fermion bound-state contributions and to extract the relevant quantity describing the 
bound states, namely the Bethe-Salpeter amplitude Xp- summary: 

— We see that a three-fermion bound state with mass M indeed gives rise to a first order pole in the total three-body 
energy P" at the bound-state energy P° cup = vP^ + M^ or, equivalently, P^ M^, P° > 0. This analytical 
dependence of Gp on P is a useful criterion to identify bound states. Note that it is just the Fourier transform of 
the function, due to the particular time-ordering (HO), which causes this singularity. 

— A further striking feature is that the Green's function becomes separable on the mass shell of the bound state, 
i.e. the dependence on the relative momenta and also in the indices for the three incoming and outgoing quarks 
separates; the product of both parts, which corresponds to the residuum of the six-point-Green's function at the 
baryon pole P° — > wp, just defines the Bethe-Salpeter amplitude and its adjoint, see also fig. ^: 



{P(..Pv\P'i^Pv) 



Xp ai 0203 {Pi,Pv) Xp {p'i^Pv)- 



(50) 



Evaluating the inhomogeneous integral equation ( pGI ) for the six-point Green's function Gp at the bound-state pole 
and using this special behavior ( pO| ) of Gp at this pole position will allow to derive the Bethe-Salpeter equation for 
the Bethe-Salpeter amplitudes and the corresponding normalization condition. This will be shown in the next section. 



-1 



lim 



G - = 



P"- M +ie / 

Baryon \ 
(bound state) 



Fig. 7. The behavior of the six-point Green's function G in the vicinity of a three-quark bound-state pole of a baryon with mass 
M: Via the adjoint Bethe-Salpeter amplitude xp as 'vertex' (right halved bubble) the three (off shell) quarks form a bound 
state (baryon), which then propagates by means of the propagator ~ (P^ — + (denoted by the threefold line) and 

finally 'decays' again via the 'vertex' given by the Bethe-Salpeter amplitude Xp (left halved bubble) into three off shell quarks. 



3 Bethe-Salpeter equation and normalization condition 

With the results of the foregoing sections, we are now in the position to derive 

— the Bethe-Salpeter equation for the Bethe-Salpeter amplitudes, which is an homogeneous integral equation de- 
scribing the bound states relativistically, 

— the normalization of the Bethe-Salpeter amplitudes. 
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This can be done simultaneously in a simple and appealing way by a Laurent expansion of the inhomogeneous integral 
equations ( |3^) and for the six-point Green's function Gp in the total energy variable around the bound state 
pole at P = P. To this end it is convenient to bring the integral equation ( ^6| ) and its adjoint (^) into the equivalent 
forms 



[Gop^ + i Kp] Gp = I, 
Gp [Gop'+iKp] = I, 



(51) 
(52) 



where the dependence on the four momentum P appears only on the left hand side. Here I is the identity for the 
operator product (35), which reads explicitly 



(53) 



and the operator Gop^ is the inverse of Gop with respect to this operator product, which thus obeys Gop^ Gq 
Go Gop^ = I. It is given by the triple product of the inverse quark propagators 



Gop'(pc,Pr,;P^,P;) = ^F + Pi + iPl) ® - Pi + ^Pri) ® S% \\P - p^) 

X {2n)U('\p^~p'^) {2n)U^'\p^-p'^). 
Equations ( ^ ) and (|5^) imply that Gp is the resolvent of a pseudo-Hamiltonian 



Gq, 



.e. 



Hp Gp — Gp Hp 



(54) 



(55) 



(56) 



In order to obtain an equation for the Bethe-Salpeter amplitudes and their normalization condition from (pq), we 
use the analytical dependence of the six-point Green's function Gp on P in the vicinity of the bound-state pole at P 
derived in the preceding subsection. Therefore we perform an expansion of the Green's function Gp and the pseudo- 
Hamiltonian Hp in the variable around the bound-state energy uip. Due to eq. ( ^ ) we find a Laurent expansion 
of the Green's function Gp beginning with the first order singularity]^. 



Gp = 



Xp Xp 



2ujp po - ujp 



d 

apo 



(P° - c^p) Gp 



O (P" - C^p) . 



and analogously for the operator Hp we have the Taylor series expansion 

d 



Hp — Hp 



QpO 



Hp 



(P°-c.p) + O (P"-c.p) 

po=wp ^ 



(57) 



(58) 



Inserting both expansions ( p7| ) and ( pq ) into eq. 
equation Hp Gp = I up to the first order: 



Hp [xpxp] {P°-Lop 
' d 



then yields the following Laurent expansion of the operator 



2cjp 
+ Hp 



QpO 

O (P° - ujp 



[(P" - ^p) Gp] 



2a;p 



9P0 



Hp 



Xp Xp 



P"=c. 



= I 



order — 1 

order 
orders > 1 
order 



(59) 



Comparing the expansion coefficients of each order in (|59|) we obtain simultaneously the equation for the amplitudes 
Xp, i.e. the Bethe-Salpeter equation, and the normalization condition. 



^ Note that with our notation P = {uJ-p, P) and P = (P", P) the Bethe-Salpeter amplitudes Xp a^nd Xp do not depend on P^\ 
as they are on shell amplitudes by definition. 
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3.1 The Bethe-Salpeter equation for three bound fermions 

The expansion coefficients in the Laurent series ( p9| ) of the order (P° — t^p)^^ yield 

Hp [XP XP] = 0. (60) 



Now the factorization property of the pole residue becomes crucial: due to the separability ( |50[ ) of the product xpXp 
the operation ( |35| ) of Hp acts only on the indices and relative momenta of xp: while xp remains unaffected, thus 
producing the Bethe-Salpeter equation for the Bethe-Salpeter amplitude xp- 

Hp XP = 0. (61) 



Here the operator product of a six-point function Hp with a three-point function xp is defined, analogous to (35), as 



[Hp Xp]aia2a3 (P^^Pv) J (Stt)^ ^ P '^^'^^as:, a'^a'^a'^iP^ i Pv' P^'Pv) XP a[a'^a'^{P^ t P,^) ■ 



(62) 



In the same fashion the corresponding Laurent expansion of the adjoint equation GpHp = I gives the adjoint 
Bethe-Salpeter equation for the adjoint amplitude^ x_p • 

XP Hp = 0. (63) 

Inserting definition ( jSq ) for Hp and multiplying by Gop, we bring the Bethe-Salpeter equation and its adjoint into 
their more conventional form 

Xp = — i Gop Kp xp, , s ,21 

with Kp = K^^> + Ky. (64) 

XP = -i XP Kp Gop, 

The three-body Bethe-Salpeter equation is a covariant eight-dimensional homogeneous integral equation in the vari- 
ables p^ = (p^, pj) and Pjj — (p5J, p^) describing the properties of bound states. It reads explicitly: 

XP (Pi^Pv) = 

S'f a,a', ilP + P« + iPv) Si C^P-Pi + kPn) S% ,3,. (iP - P,) 
^% ^(3) 



^^'^ J (2^ (2^ ^Pa[a',a',- a'^a'^a'^(Pi^ P^'^ P^Pr,) XP a'^a'^a'^^Pr,) 



J ^ ^(§Wp„) a[a',; a'^a'^iP^^P^ XP a'^a'^aM'Pv) 
+ S'p a,a', {}P +Pi + \P,) Si ,3,, (iP - P„) 

d^ ' 

^ / (2^ ^(§-P+p„2) a'^a'^; a'^ a!^ ^P^^ ^ Pi2) XP a'^a^a'^ ("5^42 + ZPv2,-p'i2 " iPm) 

+ Si a.a', ilP -Pi + iPn) Si ,,,, {\P - P,) 

/ ^(§Wp,J a^a- a^'4'(P«i'4) XP a.a'^a'^ H ^ i^m : P^l " ^^m) ■ (65) 

Recall that the two sets (p^nPrn) and {ps^2TPm) of relative momenta are related to the standard set {p^,Pri) = {p^a^Pm) 
by cyclic permutations of the quark momenta represented by the linear transformations (|30|). The Bethe-Salpeter 
equation is represented diagrammatically in fig. ^. 

^ Here the operator product is defined similar to (^) but with summation and integration over indices and momenta that 
appear on the left in Hp 
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Fig. 8. Grapliical illustration of the three-fermion Bethe-Salpeter equation for the Bethe-Salpeter amplitude xp- K^^^ and K^'^^ 
denote the irreducible three- and two-body interaction kernels, respectively. Thick arrows on quark lines indicate full quark 
propagators. 



3.2 The normalization condition 

Comparing the expansion coefficients of order — uj-p)^ in the Laurent series (Eol) gives 



Hp 



d 



gpo 



[(P°-c.p) Gp] 



2ujp 



d 



Hp 



XP Xp 



(66) 



which expresses the requirement that the product of Xp and Xp is the residue of the bound-state pole in Gp. If we 
multiply this equation from the left hand side with the adjoint amplitude xp, the first term in ( |66| ) vanishes according to 
the adjoint Bethe-Salpeter equation and we find the normalization condition for the Bethe-Salpeter amplitudes 


d 



1 XP 



QpO 



Hp 



XP 



p«= 



The full explicit expression then reads 



Xp{p'^,p'rj 



(27r)4 (27r)4 
d 



(27r)4 (27r)4 



(67) 



(68) 



XpiPi^Pn) 



2^1 



PO=u 



A priori the normalization condition provides the correct relation between the amplitudes xp ^^^d the six-point Green's 
function G. But furthermore, this additional boundary condition is essential in selecting the proper solutions xp of 
the three-fermion Bethe-Salpeter equation (|6^ ) thus providing a discrete spectrum = M'^ of bound states. 

Note that the normalization condition for the amplitudes as written in the form of eq. ( |67| ) is not manifestly 
covariant in contrast to the Bethe-Salpeter equation. But it holds in any frame since both sides of eq. ( |67|) transform 
like the time component of a four- vector, if the amplitudes transform properly under Lorentz transformations. However, 
we would like to remark here that the normalization ( |67|) may also be rewritten in explicit covariant form as 



1 Xp 



pfj. 



d 



Hp 



Xp = 2M' 



(69) 



For a diagrammatic illustration of eq. mu) see fig 



-iP 



^"-^ \ . " . nenTi 



iK' 



(2) 



+ 



\Sp®Spi8)SF 




p=p 



Fig. 9. The normalization condition for the Bethe-Salpeter amplitudes. The filled circles denote the inverse quark propagators, 
7^(2) and K^^^ are the irreducible two- and three-quark interaction kernels. 
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4 Reduction to the Salpeter equation 
4.1 Motivation and general remarks 



In principle, the Bethe-Salpeter equation (|65|) for three fermions, derived in the foregoing section, provides a suitable 
starting point for the covariant description of baryons as bound states of three quarks in the framework of QCD. Solving 
this equation for given single quark propagators Sp and interaction kernels K^^^^ and K^^\ the discrete spectrum of 
states is then determined by the normalization condition (^9|). However, an exact solution of the Bethe-Salpeter 
equation within the framework of QCD is impossible, since the quark propagator Sp and the irreducible interaction 
kernels if*^^^ and K^^^ are only formally defined in perturbation theory as an infinite sum of Feynman diagrams. 
Moreover it is unclear, which particular approximation will provide quark confinement in hadrons. 

But even if the exact kernels and propagators were known in QCD, the dependence on the relative energy (or the 
corresponding relative time) variables leads to a complicated analytic pole structure, which so far could be treated 
rigorously only in the case of two scalar particles interacting through a (massless) scalar exchange (the so-called 
Wick-Cutkosky model, see Thus, the use of general two-quark and three-quark interaction kernels, that 

depend on the relative energy variables, leads to serious conceptional and practical problems. To our knowledge the 
only attempt to solve (approximately) a full four-dimensional three-quark Bethe-Salpeter equation in Euclidean space 
has been performed by Meyer and Bohm and subsequently by Kielanowski ||l5| and Falkensteiner |l6j in an 

approach where baryons were considered as extremely strongly bound systems of three quite heavy constituent quarks 
m ^ IGeV, that interact via harmonic oscillator interactions, so that a solution can be obtained by an expansion in 
powers of ^. However, from a modern point of view, the crude approximations and especially the large constituent 
quark masses are questionable and not suited for phenomenologically successful applications. 

Thus, the use of the full eight-dimensional Bethe-Salpeter equation is of rather limited practical value and the lack 
of a confinement kernel that could be rigorously derived from QCD anyhow requires an appropriate phenomenological 
parameterization: so far, the only ansatz that can give a realistic description of the quark confinement and thus can 
account for the gross features of the whole baryon spectrum up to highest orbital excitations, is the nonrelativistic 
quark model, which uses static two- and three-quark potentials. 

For these reasons we will not treat the full three-quark Bethe-Salpeter equation. Instead we try to eliminate the 
difficult relative energy dependence in order to get a six-dimensional reduction of the full eight-dimensional Bethe- 
Salpeter equation, the so-called Salpeter equation IT^, with the aim to obtain a framework that is still covariant. 
At the same time we want to keep as close as possible to the quite successful nonrelativistic quark potential model 
in order to obtain at least this model as a non-relativistic limit. In this spirit, a covariant quark model for mesons 
based on the instantaneous gg-Bethe-Salpeter equation has been developed already and has been successfully applied 
to the calculation of mass spectra and various transition matrix elements up to high momentum transfers, see 
|l9|] . To extend this model for calculations of baryons, we make the same simplifying assumptions and approximations 



in the three-quark Bethe-Salpeter equation (sect. 4.2): The full quark propagators Sp are assumed to be given by 
their free forms with effective constituent quark masses. Moreover, the kernels if^^-* and K'-^'> are approximated by 
effective interactions that are instantaneous in the rest frame of the bound state, which thus corresponds to the 
neglect of retardation effects. We should mention here that the instantaneous approximation can be formulated in a 
frame independent way psf , so that formal covariance is preserved, which becomes important for the calculation of 
transitions between baryon states, where at least one of the baryon has to be boosted. 

In the meson case these approximations allow for a direct and straightforward reduction to the gg-Salpeter equation 
by an analytical integration over the relative energy variable, since the connected instantaneous gg-kernel 
cuts the whole relative energy dependence of the Bethe-Salpeter equation. The same applies also to the three-quark 
Bethe-Salpeter equation, if only an instantaneous, connected three-quark kernel K'--^") is taken into account and two- 
particle kernels are neg lected (/i:(2) ^ 0). In this case the Salpeter equation can be formulated in a concise Hamiltonian 
form with some characteristic projector properties that reduce the number of independent functions necessary to de- 
scribe a baryon state. For the sake of conceptual simplicity such an approach has been used in our former investigations 
pq,E^Pq], where all kinds of interactions have been parameterized in a kind of collective instantaneous three-body 



kernel. In section 43 we will first give a summary of the reduction procedure in this simple and instructive case and 
discuss the specific structure of the resulting Salpeter equation. 

However, as soon as genuine two-quark kernels K^-^') are considered, new difficulties arise since the two-body terms 
are unconnected within the three-quark system: despite an instantaneous approximation of K^^') there remains a 
relative energy dependence due to retardation effects of the third non-interacting spectator quark, which is off-shell in 
general. In this respect the elimination of the relative energies is technically and conceptually much more involved and 



an enhanced reduction procedure is needed. In section 4.4 we give a procedure that nevertheless allows for the reduction 
to a Salpeter equation. The crucial point is the existence of a genuine instantaneous connected part of the interaction 
K'^^\ right fr om start. In our model this part will be given by a convenient form of a static three-body confinement 
potential that must be present for all baryon states in all sectors due to the confinement hypothesis. Recasting the 
Bethe-Salpeter equation into a more convenient form with all two-particle effects collected into a six-point Green's 
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function thus provides a similar reduction procedure as in the case of vanishing two-body interactions. Extending a 
kind of quasi-potential approach as it was first proposed by Logunov and Tavkhelidze ||29|| for the equal-time Green's 
function of two scalar particles, all effects of the unconnected two-body interactions can then be transformed into an 
effective instantaneous potential that adds to the genuine three-body kernel K'-^^ and we finally end up with a reduced 
equation that exhibits the same expedient projector structure as in the case where the dynamics of the quarks is given 
by an instantaneous three-body kernel alone. The effective potential, however, consists of an infinite perturbation series 
of time-ordered Feynman diagrams, which needs to be truncated for explicit calculations. In the subsequent sect. |^ 
we will analyze the structure of the resulting baryon Salpeter equation and its corresponding Salpeter amplitudes in 
detail: a remarkable substantial property of our covariant Salpeter approach will turn out to be that it exhibits a 
one-to-one correspondence with the states of the nonrelativistic quark model. 



4.2 Approximations 

In order to transform the Bethe-Salpeter equation into a solvable integral equation several simplifying approximations 
have to be made. To start, we follow the prescription of and assume free quark propagators and instantaneous 
interaction kernels. 



4.2.1 Free quark propagators 

First, we make the assumption that the full quark propagators can be approximated by the usual free fermion propa- 
gators with effective constituent quark masses for each quark|^ 

pi - mi + le 

This approximation is consistent with the picture of a hadron mainly built out of constituent quarks analogous to the 
non-relativistic quark model. The effective constituent quark masses rui enter as free parameters in our model. 



4.2.2 Instantaneous approximation 



Moreover, we choose the irreducible two- and three-body interaction kernels to be instantaneous in the rest frame of 
the baryon, meaning that in the center-of-mass system there is no dependence on the relative energy variables and 



P=(M,0) 



(71) 



(2) 



P={M,0) 



(72) 



This approximation corresponds to the neglect of retardation effects in the rest-frame. To preserve the formal covariance 
of the Bethe-Salpeter equation, however, we need a covariant description of the instantaneous approximation which 
holds in any arbitrary reference frame of the bound state. Wc follow an idea of Wallace and Mandelzweig j25j and 
introduce for arbitrary time-like total four-momenta P, P^ > 0, the following covariant decomposition of any four- 
dimensional four-vector p, 

P 

P = P\\ —F^+Pi- (73) 



VP2 

into components parallel and perpendicular to the total four-momentum P 

{P,P) {P,P\ 



P\\ 



rp2 



p±:^p 



p2 



P. 



(74) 



This is a decomposition into a time-like vector p|| P/V P^ and a space-like vector p± which effectively is three- 
dimensional in content. Now the instantaneous approximation, which has been formulated in eqs. ( [7l| ) and ( [72| ) within 

^ For a simplified notation we suppress the explicit flavor- and color dependencies for the moment. 
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the center-of-mass frame of the three-body system, can be formulated in any reference frame (which is specified by the 
four- momentum P): we assume that the kernels do not depend on the time- like parallel components of the relative 
momenta, i.e. P5j|,p,, ,p' , but only on the space- like perpendicular components: 



Kp\p^,p,{, P^,P^) = V^^\p^^,p^^; p^^,pIj^), 



(75) 



K 



(2) 



For interaction kernels of this type we have 



d 



P'^,Kf{p,^p,- p'^,p'^) 



P=P 



= 0, 



(76) 



(77) 



p=p 



P=P 



P=P 



and, consequently, these give no direct contributions to the normalization condition 
plitudes. In the rest frame of the baryon where P = P = {AI, 0) we find 



P\\ =P 



and 



P± = (0,P) 



(78) 



f) for the Bethe-Salpeter am- 



(79) 



so that the covariant formulation of the instantaneous approximation given in (|7^) and ( [tg] ) indeed recovers the con- 
ditions (^Tj) and ( |7^ ) in the center-of-mass frame. 

In the two-fermion case [T^ ] it was shown that the assumptions of free quark propagators and instantaneous 
interaction kernels are sufficient to completely eliminate the dependence on the relative energy dependence in order 
to arrive at a reduced equation which can be solved with standard techniques. In the three-fermion problem, however, 
this is in general not possible, unless we consider systems without two-body interactions. In the more general case 
new difficulties arise from the property of the two-body terms that these are unconnected within the three-body 

(2) 

system. Despite the instantaneous approximation of the two-particle interactions, the kernel ifp=(Mo) defined by eq. 
( p^ ) remains retarded, since (in the CMS) it maintains the dependence on the relative energies p° and p° due to 

(2) 

retardation effects of the third noninteracting spectator quark which is off-shell in general. Accordingly, in Kp^j^jyj q) 
this remaining relative energy dependence is given explicitly by the inverse single quark propagator of the spectator 
together with its four-momentum conserving (5-function: 



Kp\p^,p,j; p'^.p'rj) 



P={M,0) 



of(123) 



Pm) 



P={M,0) 



(80) 



Thus, the consideration of unconnected two-particle terms in the three-body Bethe-Salpeter equation makes a reduction 
technically much more involved, despite the instantaneous approximation of the two-body kernels. With regard to the 
goal of finding a convenient reduction procedure it is therefore instructive to consider first the conceptually much easier 
case of vanishing two-particle kernels, where the dynamics of the quarks is determined by a connected instantaneous 
three-body kernel alone. In this case the reduction of the eight-dimensional Bethe-Salpeter equation to an equivalent 
six-dimensional equation - the so-called Salp eter equation - is straightforward (as in the two-fermion case with a 
connected instantaneous two-body kernel p8[ ). 



4.3 The reduction without two-particle kernels 

Neglecting the irreducible two-particle interaction kernels, i.e. K^^'' — 0, and taking only an instantaneous three-body 
kernel (|75[ ) into account, the Bethe-Salpeter equation and its adjoint in the center-of-mass frameQ of the baryon with 

* Due to the formally covariant formulation (^^ of the instantaneous approximation of the irreducible three-body kernel 
(which preserves the formal covariance of the Bethe-Salpeter equation), it is sufficient (and convenient) to go into the center- 
of-mass (CMS) frame 
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P = (M, 0) = M are given by 



XM = -i GoM V^^'' XM, 



(81) 
(82) 



The crucial point is now that V^^^ being instantaneous truncates the p^, p'^ dependences of the Bethe-Salpeter equa- 
tions (p^) and (p2|). This has the following consequences: 



1.) The p^, p^ integration within the operator product V'^^^ xm on the right hand side of eq. (81) acts on xm 
directly and thus can be used to reduce this eight-dimensional Bethe-Salpeter amplitude to a six-dimensional amplitude 
<?M, *-e- in detail 



XM 





(27r)4 






(27r)3 


(27r)3 


d'p'^ 





V''^\p^,Prj; P'^, P'r,) XM{p'^,p'r,) 



(2^ (2^ ^ (P^'P"' P^'P") J 



XM{p'e,p'n) 



(27r)3 (27r)3 



1/(3) 



M 



(P«,Pr,)- 



(83) 



Consequently, there remains a six-dimensional integral operation^ of y(3) on the reduced six-dimensional amplitude 
^Af, which is the so-called Salpeter amplitude: 



^Af(pe,p,,) 



-{^ ^ XM {{plp,l{P>n)) 



(84) 



In the same way one proceeds with the operator product Xm ^'^^^ in the adjoint Bethe-Salpeter equation (^), i. 



Xm V^'^ 



(P«,P.)): 



which accordingly defines the adjoint Salpeter amplitude 

dpi dp"^ 



<^A/(PC,P»?) 



Xm ((P?,P?),(P°,P»j)) 



(85) 



(86) 



(87) 



27r 2tt 

2. ) Inserting (|8^) and ( ^5| ) into the Bethe-Salpeter equations (^ij) and (^^, respectively, we have 

XM = -i Gom v'-^^ <?m, 
Xm = -i^M V^^^ Gom, 

which gives a prescription how to reconstruct the full Bethe-Salpeter amplitudes from the Salpeter amplitudes for 
any on-shell total momentum. Consequently, in the instantaneous approximation it is sufficient to know the reduced 
six-dimensional Salpeter amplitudes <?Af and <?m to get the full eight-dimensional Bethe-Salpeter amplitudes xm and 
Xm^ ^-S- the solutions of the Bethe-Salpeter equation (^ and (^), respectively. The next step is to get an equation 
which determines 'Pm and (I'm- 

3. ) As shown in eqs. ( ^ ) and (|85|), the quantities 

rM{p^,P-n) [Go a/ Xm] iPi,Pv) = [V^^^ xm] {p^,Pn) = [^'^^ ^m] (P5,P.,) = A.f(P5,P,,), ^^^^ 

rM{Pi,Pv) := [xm GqIj] (Pi,P,,) ^ [xm y''^^] (Pi^Pv) = [<^M V^^'^] {Pi,Pr,) = TM{Pi,Pn), 

^ Notice that we do not introduce a new product notation for this six-dimensional integral operation to distinguish it from 
the eight-dimensional one. The difference between the two products should be obvious from the context in which they are used. 
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which are usuaUy caUed amputated Bethe-Salpeter ampUtudes or three-quark vertex functions, do not depend on 
the relative energies p^'^ and p,,'^ in the center-of-mass frame of the baryon. Consequently, the analytical dependence of 
the Bethe-Salpeter amplitudes xm — Com and Xm — Gq jif on the variables and Pri^ stems exclusively from 
the triple tensor product Gqm of the free quark propagators. This enables us to reduce the eight-dimensional Bethe- 
Salpeter equations for the Bethe-Salpeter amplitudes to six-dimensional integral equations for the Salpeter amplitudes 
by integrating out the Pi,*^ dependence on both sides of eqs. ( ^7| ) and (^8|). The Bethe-Salpeter amplitudes on 
the left hand side reduce to the corresponding Salpeter amplitudes and on the right hand side the relative energy 
integration affects merely the free propagator Gomj '-G- detail 



^M(PC,Pr,) 



"^4 '^P° (( ^ / ^^ 

^ ^XM((p^,Pe),(p^,Pr,)) 



/ J2^ (GoM>(Pe,P,; Pe^P^) [^^'^ ^m] (p^,p,/) 



(Gom) V^^'^<?Af (P«,P^) 



(90) 



Thus, we end up with the so-called Salpeter equation and its adjoint for the Salpeter amplitudes (I>m and <?a/ 

= -i (Gom) V^^^ <Pm, (91) 
= -i^M V^^'^ (Gom)- (92) 

Here we introduced the notation 



(A)(p„p,; p^,p;) ■■=J'4'4 J'4'-i 



(93) 



for the six-dimensional reduction of any eight-dimensional six-point function A. Accordingly, (Gom) = {SpiS)Sp^Sp) 
denotes the reduction of the free three-quark propagator Gom defined in eq. (|26|). Due to the approximative choice of 
bare quark propagators with effective constituent quark masses, the analytical structure of Gqm ™ the relative energy 
variables p'^ and pjj is rather simple and consequently, the p^, p° integration in (Go a/) can be performed analytically 
by applying Cauchy's theorem. To this end it is convenient to use the followi ng p artial fraction decomposition of the 
free one-particle propagators into positive and negative energy contributions |30[, 



S'pip^) = 



cji(pi) + ie 



-tjj(pi) - ie 



which isolates the pole positions in the energy variable located at the relativistic on-shell kinetic energies 

t^i(Pi) := 



(94) 



(95) 



of the quarks. The operators yl^(pi) are the projectors onto positive and negative energy solutions of the free Dirac 
equation, written explicitly as 

a;,(Pi) i±H,ipi) 



where Hi is the usual free single particle Dirac-Hamiltonian given by 

Hi{pi) := 7° (7-pi + rrii) = cc • pi -|- /3. 



(96) 



(97) 
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Performing the pJJ integration we obtain the three- fermion Salpeter propagator: 

(Gom)(p«,p»,; p^,p^,) = 



At{pi)(g>A+{p2)(g>A+{p3) 



ipi)(g)A^ (P2)®^3 (P3) 



M - cji(pi) - tJ2(P2) - ^sCps) + ie M + cji(pi) + cj2(p2) + t^3(P3) - ie 



X 7 



.7° 7° (27r)3^(3)(p^_p,) (2^)3^(3)(p^^_p/) 



(98) 



with Pi — Pi(p{, Pr,) defined as in eq. ( |2C| ) with P = Pi +P2 + P3 = 0. Notice the remarkable property that due to the 
pole structure of Gqm in the relative energy variables and p^, the residue theorem merely provides the projectors 
onto purely positive-energy and purely negative-energy three-quark states. All mixed components vanish! 
Finally, the Salpeter equation (^Tj) in the case of vanishing two-quark kernels reads explicitly 



^M(PC,Pr, 



/l+(pi)®yl+(p2)®ylJ(P3) 



A^ (Pl)(g)^2 (P2)0^3 (P3) 



M - t^i(pi)-tJ2(p2)-t^3(P3) + ie M + cji(pi)-|-cj2(p2)+t^3(P3) - ie 



7 (X) 7 



7 



Pi 



(27r)3 {2t:Y 



y(3)(. 



M 



(Pc,p;)- 



(99) 



Thus, we have seen that in the case, where the dynamics of the three quarks (fcrmions) is described by an instan- 
taneous, connected three-body kernel alone, the reduction of the full eight-dimensional three-fermion Bethe-Salpeter 
equation to the six-dimensional Salpeter equation (in the CMS) is straightforward. The Salpeter equation is equivalent 
to the full Bethe-Salpeter equation since eq. (|8^) allows an exact reconstruction of the Bethe-Salpeter amplitude xm 
from the solution (I'm of the Salpeter equation in the rest frame. Finally, the formally covariant framework provides the 
possibility to obtain the amplitude xp in frame with = by a kinematical Lorentz boost of the rest-frame 
amplitude xm- 



According to eq. (^9|) we find the remarkable fact that the reduction in the case of pure instantaneous three- 
body kernel leads to certain projection properties for the Salpeter amplitudes which effectively reduce the number of 
independent functions necessary to describe a baryon state. Let us continue our discussion with some investigations 
of this particular structure of the Salpeter equation . 



4.3.1 The projector structure of the Salpeter equation 

Due to the energy projectors appearing in the Salpeter propagator (Gqm), the Salpeter amplitudes are eigenstates of 



the Salpeter projectors 

yl(p^, p^) := yl^ (pi) ® A+{p2) ® ^^(P3) + ^r(Pi) ® (P2) ® (P3), (100) 

^(P?, P.)) - 7° 7° 7° ^(P«, P„) 7° ® 7° 7°- (101) 
which project onto the subspace of purely positive and negative energy components, i.e. 

$M = A $M = + <Pm'-, (102) 

= ^mA = + ^m"" . (103) 

and accordingly the Salpeter equation only involves the amplitudes 

$+++ := yl+ (g) yl+ ® A+ <?M and A' <g) A' iSi A' ^m, (104) 



whereas all mixed components such as ^^J/" vanish. This property reduces the Salpeter amplitudes effectively to only 
16-component functions of the six variables p^, Ptj, in contrast to the full (in Dirac space) 64-component Bethe-Salpeter 
amplitudes, which are functions of eight variables. This projector structure implies that for the dynamics of the three 
quarks in the bound state (baryon) not the full structure of the instantaneous three-body kernel V^^"^ is relevant, but 
only its projected part 

Vi'^(p?,p^; Pe,p;) :=:4(p5,p„) F(3Hp^,p^; p'^.p'^) Aip'^,p'^). (105) 
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Consequently the residual part :— y*-^-* — , which describes the coupling to the mixed energy states, plays 
no role for spectroscopy {i.e. the determination of bound state masses), although they become relevant for the re- 
construction of the full Bethe-Salpeter amplitude xm according to eq. ( p7| ) and thus can contribute when calculating 
various transition matrix elements [pll. 




Fig. 10. Time ordered graphs of an instantaneous three-body interaction which contribute to the three-quark propagation in 
the Salpeter equation. The instantaneous three-body kernel is represented by the dashed line. 

In the language of time-ordered perturbation theory this means that the instantaneity of the kernel prevents the 
inclusion of single and double Z-graphs which correspond to the mixed components ~ , A ^, . . ., etc. of the 

interaction kernel. However, compared to a nonrelativistic ansatz, where all three quarks propagate forward in time 
(corresponding here to the components ^ A^'^^), the Salpeter equation takes into account also those diagrams, where 
all three quarks propagate backwards in time (triple Z-graphs corresponding to the components ~ A and their 

coupling to components ~ via A V^^'' /!+"'""'"), as shown in fig. We want to remark here that the 

appearance of these negative energy components in the Salpeter equation is connected with the particle- antiparticle 
symmetry due to the CVT invariance of the underlying relativistic field theoretical framework. We w ill come back 
to this characteristic feature of the Salpeter equation and discuss it in some more detail in sect. 5.2 after we have 
taken also the two-particle interactions into account. The importance of the negative energy contributions depends on 
the energy denominators M =p {uji -\- u>2 + (j-i^ — ie) of the positive and negative energy components in (99) as can be 
illustrated by the following two extreme cases: 

— For small binding energies, i.e. M Ki mi + m2 + and |pi|/mi <C 1 one has 



1 



1 



M + LiJ\ + UJ2 + U}^ M — LOl — U}2 — OJ^ 



(106) 



such that the negative energy component in ( |99| ) becomes negligible compared to the positive component and one 
is led to the so-called reduced Salpeter equation. 

For deeply bound states, i.e. M <^ mi + m2 + ms, both components are of equal order of magnitude: 



1 



1 



M + uji + LU2 + ^^3 M ~ uji — UJ2 ~ toy, 



(107) 



In our case of baryons as a bound three-quark system we should definitely be rather far away from the limit of deeply 
bounds states. Nevertheless the negative energy term of the Salpeter amplitude might contribute to a certain amount. 



4.3.2 Hamiltonian formulation of the Salpeter equation 

The special projector structure in connection with the particular energy denominators M =p (oji + 0)2 + 0J3 — ie) allows 
for the formulation of the Salpeter equation in Hamiltonian form, i.e. as an eigenvalue problem 



M- 



H <Pm = M with A(pR 

Here we define the Salpeter Hamiltonian H by 

[n<pM] (P?,P,,) Wo(PC,Pr,) '^M(PC,Pr,) 

+ [ylf(pi) «) A+{p2) (E) A+{p3) + /Ir(Pi) ® (P2) <E) A^{p3) 



(108) 
(109) 



X 7 (g) 7 
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where the free three-fermion Hamiltonian 

Ho(P5,Pr,) i/i(pi)®I®I 
represents the relativistic kinetic energy operator. 



I(8)i?2(P2)' 



1(g) 1(g) i73(p3) 



(110) 



Of course, a similar representation of the adjoint Salpeter equation, which determines the adjoint amplitude ^m, 
can also be found. Note however, that both equations are not independent, but even are equivalent, since there is a 
general]^ interconnection between the Salpeter amplitude and its adjoint 'Pm, which in momentum space reads: 



<pM{Pi,Pv) = -'^L(P?,Pr,)7"®7''®7"- (111) 
To be consistent, one has to require: If is a solution of the Salpeter equation (|108[), then <1>m, as defined by 



relation (111), has to be a solution of the corresponding adjoint Salpeter equation (and vice versa). Using Hq = Ho 



and A] = Ai, one easily shows that this equivalence of the Salpeter equation (lOS) and its adjoint implies the following 
condition for the interaction kernel y'^' in the CMS: 



i^^'^Pop;; p?,p.) 



it 



7 (g 7 



(112) 



4.3.3 Normalization of Salpeter amplitudes - Scalar product 



The normalization condition ( |67[ ) of the Bethe-Salpeter amplitudes, which reads in the center-of-mass frame with 
P = (M, 0) = M 



Xm 



QpO 



(113) 



P"=M 



induces a normalization condition of the corresponding Salpeter amplitudes . The instantaneous three-body kernel 
V^'^^ has no explicit energy dependence and thus gives no contribution to the norm. Using the representation xm = 
Gqm rM and Xj^^ = Fm Gqm of the Bethe-Salpeter a mpli tudes, where the vertex functions Fm and Fm defined in 
( p9| ) do not depend on the relative energies p", eq. (113) becomes 



2M = -i Fm Go 



AI 



d_ 

dM 



T77T Gom 



GoM A/ — F]^[ ( — i Go 



3M 



d_ 

dM 



^77 ^Om 



Go 



OM > rM 



(114) 



Here the angled brackets (. . .) indicate the internal integration over and pJJ which is used to reduce the enclosed 

operator according to eq. (p3|). With Gom {d/dM Gom) Gqm — —d/dM Gom this reduced operator may be 
rewritten as the derivative of the Salpeter propagator (p8|) and we obtain 



-i Go 



M 



d 
'dM 



Gqm 



OM 



'A[ ^^'''^ ^ 7° ® 7° ® 7° (Gom). 



(115) 



Substitution into eq. (114) and replacing the vertex functions according to the relations (Go7\/) — 'Pm and 

AI 



Fm {Gom) = '^ai ~ ^^1/ 7" ® 7° ® 7° then yields the following normalization condition of the Salpeter 



amplitudes <Z> 



M- 



(116) 



Thus, the solutions <?m of the Salpeter equation have to be normalized according to the usual £^-norm just like the 
solutions of the ordinary nonrelativistic Schrodinger equation. This norm induces a positive definite scalar product 
for arbitrary amplitudes ^1 and ^2 that are restricted to positive and negative energy components, i.e. <Pi — A<Pi: 



/d^ d*^ 
H '^*la,a,a,{Pi,Pn)'^2a,a,as{Pi,Pn) 



(117) 



^ i.e. the interconnection (111) between $m and $m is independent of the so far considered assumption of vanishing two-body 
kernels and other approximations of the Bethe-Salpeter equation. 
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Hence, the normalization condition (f f6) for solutions of the Salpeter equation is then given as 



(118) 



We want to remark here that a similar treatment of the fermion-antifermion (or the two-fermion) system does not lead 
to a positive definite scalar product, owing to a relative sign between the positive and negative energy contributions, 
see il8i. 



Note that the Salpeter Hamiltonian H is hermitean with respect to the scalar product (117), i.e. 

{<Pl\H<p2) = {H<Pl\'p2) 



(119) 



which is a direct consequence of the condition (112) on V'^"^^ resulting from the interconnection (111) between <1>m and 
<1>M- This guarantees, as in the case of the ordinary Schrodinger equation, two important consequences, namely: 

— The eigenvalues M of H, i.e. the three-fermion bound-state masses are real, i.e. M* = M. 

— Salpeter amplitudes <?Mi and <?m2 belonging to different eigenvalues Mi M2 are mutually orthogonal, i.e. 



4.4 The reduction with genuine two-particle kernels 

Now let us come back to the general case we are in fact interested in, where in addition to the instantaneous three- 
body kernel V^^\ the dynamics of the quarks is also determined by the unconnected instantaneously approximated 
two body-terms given by ( [72[ ) and ( |80| ) . Referring again to the formal covariance of the instantaneous approximation 
as before, we choose for these considerations the three-body rest-frame with P = (M, 0) = M. The Bethe-Salpeter 
equation and its adjoint then read 

XM = -i Com V^^^ xm ~ i Gqm Xm. (120) 



Xm 



4 Xm V^^^ Gom - i Xm K^m Gom, (121) 



and now the difficulty stems from the circumstance that due to the second term on the right hand side of (120) 



(2) 

which contains Kl,j , the relative energy dependence in the Bethe-Salpeter equation can no longer be separated 
and thus, from the outset, the reduction cannot be performed. Nevertheless, we still can take advantage of the 
fact that the p^, p^ dependence at least is cut by the first term, due to the instantaneity of V^^\ Recasting the 
Bethe-Salpeter equation into a more convenient form, this feature will in fact provide a possibility to perform a 
reduction, as we will see in the following discussion. But let us emphasize that the way of how to perform the 
reduction and consequently the final form of the Salpeter equation is not unique, although the various resulting 
reduced equations are formally equivalent. However, in practice, even the reduced equations are not solvable in general 
so that further approximations are indispensable and thus the different reduced equations become practically non- 
equivalent. Therefore, the reduced equation in its full exact form should, right from start, have an expedient canonical 
structure allowing further approximations to be made in a systematical way. Referring to this we will orientate our 
considerations according to the instructive c anon ical fo rm o f the Salpeter equation as given in the previously discussed 



case of vanishing two-body kern els b y eqs. (108) and (|l09|). Before we present this specific method for the reduction 



of the Bethe-Salpeter equation ( 120 ) and its adjoint ( |121D in practice, let us generally discuss in a first attempt. 



how in principle it becomes possible to reduce the eight-dimensional three-fermion Bethe-Salpeter equation to an 
equivalent six-dimensional Salpeter equation, provided that the full interaction kernel Kp contains at least one 
connected instantaneous part, as given in our case by the instantaneous three-body kernel, 

what changes at all in the structure and the properties of the Salpeter equation and thus of the Salpeter amplitudes 
<?M in comparison to the case discussed previously, where the dynamics was given by an instantaneous three-body 
part alone. 



4.4.1 A first attempt - concepts, ideas and problems 

In a first attempt we now want to sketch a procedure showing that a reduction of the Bethe-Salpeter equation ( |120D can 

indeed be achieved, utilizing that y*-^-* cuts the relative energy dependence in one term of the Bethe-Salpeter equation. 

— (2) 

The crucial idea and concept of this procedure is to get rid of the problematical second term — i Gqm K j^j xm 
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appearing on the right hand side of eq. (|120| ), where the unconnected two-body term K j^,f a cts on the Bethe-Salpeter 



ampHtude xm directly. This can be reached by recasting the Bcthe-Salpeter equation (12C) in the following manner. 
First we separate the terms of the Bethe-Salpeter equation into p^, dependent and independent parts, as follows: 



1 , .77(2) 



XM 



(122) 



Remember that — i V^^^ xm on the right hand side indeed has no relative energy dependence due to eq. (^9|). Now let 



US introduce the resolvent t/j^^ of the operator [Gq"^^^ 



(2) 

i K ] appearing on the left hand side of the eqs. (122) 



I.e. 



77(2) 

yM 



Gq ^ 



M 



•1?(2) 



Gq ^ 



M 



• 1^(2)1 -(2 
1 A 



M 



<M 



(123) 



_(2) 

This Green's function Q j^j is the solution of the inhomogeneous eight-dimensional integral equation 



7^(2) 



T^C^) 7^(2) 



Gm — Go A/ — i GoM i^jvf Sm 



(124) 



and thus describes, apart from the free propagation Gqm, also the propagation of the three quarks via the unconnected 

two-particle interactions alone. Multiplying eq. ( |122| ) by this resolvent Q^^^ we obtain the Bethe-Salpeter equation in 

a form similar to (|8l|), i.e. the case where we neglected the two-particle forces, but with Go a/ now replaced by G^^i 
which additionally collects all remaining retardation effects concerning the unconnected two-quark interactions within 
the baryon: 

Xm^-iGm V^'^ XM. (125) 

This form enables us again to exploit the crucial property of the instantaneous kernel V^^'^ to separate the dependence 
on the relative ener gy v ariables p^'^ and Consequently we can proceed in the same way to reduce the Bethe- 
Salpeter equations ( 125 ) as we did when reducing eq. (^) in the case of vanishing two-quark kernels. According to 
eq. ( p3| ) the eight-dimensional integral operation y'^^^ on xm on the right hand side of eq. (125) can be reduced to a 
six-dimensional operation on the Salpeter amplitude <Pm, 



^''^ xm] {Pi,Pn) = k^'^ <^Af] (Pe,P„)- 



(126) 



This (in principle) provides us again the possibility to reconstruct the full eight-dimensional Bethe-Salpeter amplitude 
Xm from the Salpeter amplitude <Pm according to 



XM 



(127) 



and shows that the analytical p^, dependence of the Bethe-Salpeter amplitude is completely determined by the 



analytical structure of Gm these variables. Thus, performing the integration on both sides of eq. (127), the 



(2) 

Bethe-Salpeter amplitude xm on the left reduces to the Salpeter amplitude and on the right hand side only G m 

(2) 

is affected by this integration and reduces to {GJ^.f ). Consequently, we finally end up with the reduced equation which 
determines the Salpeter amplitude <Pm, ?-e. 



M 



-1 (G'i}) V^'^ <Pm. 



(128) 



All the difhculties, arising from retardation effects due to the unconnected two-body terms, are now transferred to the 

(2) (2) 1 I 

reduction {G m) of Gm ■ Corresponding to the inhomogeneous integral equation (124) this reduction of the Green's 

function is determined by 



7^(2) 

Gm 



= {Gqm) - i {Gqm G 



T7(2) 7^(2) 



'M 



= (G, 



OM) 



Gq 



M 



G, 



OA/ 



Go 



M 



-iK 



G, 



OAf 



-i K 



G, 



OM 



(129) 



Thus, we have shown that, even with unconnected two-particle kernels, it is in principle possible to reduce the 
eight-dimensional threc-fcrmion Bethe-Salpeter equation to a six-dimensional Salpeter equation, provided we choose 
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at least one part of the full interaction kernel to be instantaneous. Due to the interconnection (127) of the full eight- 
dimensional Bethe-Salpeter amplitude xm and its six-dimensional reduction, i.e. the Salpeter amplitude 'Pm, the 



Salpeter equation is equivalent to the full Bethe-Salpeter equation since eq. (127) provides in principle an equally 
exact reconstruction of the Bethe-Salpeter amplitude. Unfortunately, the analytical structure of the Green's function 

Qj^j in the complex planes of the relative energy variables and is rather complicated, so that in practice neither 



^(2) 



its reduction (Gm) required for solving the Salpeter equation (128), nor Qj^j itself, required for the reconstruction 



^ ^ ^ ^ ^ ^ (2) 

(127) of the Bethe-Salpeter amplitude, is manageable in its full, exact form. The determination of {G m ) for example 

(2) 

requires in principle the calculation of an infinite number of reduced diagrams due to the Neumann series of Q j^j , 

(2), 



see 



eq. ( |129| ). We do not want to bother about that complexity at the moment and first consider {G j^j , 
up to the Born term, 



M 



-i (Go 



Gqm aIj g, 



M 



in eq. (128) only 
(130) 



in order to inspect what changes basically in the structure of the Salpeter equation and the corresponding Salpeter 
amplitudes. A tedious but straightforward calculation, using the residue theorem for performing the relative energy 
integration, yields the following contributions to the reduced Born term: 
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For the sake of clarity we suppressed partially the explicit coordinate dependences by using the more compact notation 

Af = yl^(pi), LUi = Wi(pi) and Af = Af{p[), uo'^ = tJi(pJ). The time-ordered Feynman graphs corresponding to the 

six different terms in eq. ( |131[ ) are shown in fig. |ri|. 

The first two terms (a) and (b) have the same projector structure and corresponding energy denominators as the 

Salpeter propagator (Gom) a'^d thus are of a similar form as the reduction of a genuine instantaneous three-body 

interaction diagram. The decisive difference to the previously discussed case is due to the occurrence of the mixed 

energy components (-1- H — ), ( h), etc. in the remaining four terms (c) - (f), which result from retardation effects 

(2) 

of the unconnected two-particle interactions. In other words: The propagator {Q ) which has been substituted for 
(Gqm) (compared to the case of neglected two-body kernels) does not exhibit the particular projection properties of 
(Gqm) I the restriction to purely positive and purely negative energy components only. This implies (compare the 
discussion in subsect. 4.3.1): 



— The Salpeter amplitudes '^m are no longer eigenstates of the Salpeter projector A, but also possess mixed energy 
components according to the terms (e) and (f) in (131). 

— In connection with the unconnected, retarded two-particle kernels, also the 'residual' part vf'^ = F^^) _ yf'^ of the 
instantaneous three-body kernel '-3) ^ that couples to the mixed components, now contributes to the threc-fcrmion 
bound state and therefore to the spectroscopic results. 
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Fig. 11. Time ordered grapiis of tlie reduced Green's function {Qm ) up to the Born term (only the interactions in the quark 
pair (12) are shown.) The instantaneous two-body interaction between two quarks is repre sente d by the vertical dashed lines 
and dots on the affected quark lines. The Born graphs correspond to the expressions of eq. (131). 



Note however, that (assuming weakly bound states) the important dominant terms of this Born contribution are 
given by the purely positive and negative contributions (a) and (b) due to their particular structure of the energy 
denominators: 

— In the case M > the terms (c) and (e) are suppressed with respect to the dominant term (a), since {M — lo\ — 
(jji — 3> (wi + Cci2 + ti^i + tjj) All other contributions have denominators (A/ + + + "-^s) and thus are 
anyhow suppressed. 

— In the caseH M < the dominant term is (b), whereas (d) and (f) are suppressed against (b), since (M + cji + 
102 + ^s)"^ 3> + clI2 + ^1 + "^2)^^ All other contributions have denominators (Af — cji — ti>2 ~ "^a) and hence are 
anyway negUgible. 



The more complex structure of the Salpeter equation ( |128| ) restricts its appHcabiUty to explicit three-body bound- 
state calculations: Due to the explicit appearance of the additional mixed energy components the formulation as an 

eigenvalue problem in Hamiltonian form such as in the case of a pure instantaneous three-body kernel is no longer 

(2) 

possible. Moreover further approximations of the reduced Green's function (C/jy,^ ) are indispensable, which gives rise to 

(2) (2) 

the question of how to approximate (^Q m ) systematically. One would expect that a perturbative approximation of m ) 

(2) 

simply by cutting the Neumann series of Q at finite order {e.g. in the so far considered Born approximation), would 

not be sufficient to describe accurately the effects of the two-particle interaction within a three-body bound state, as 

(2) 

e.g. two particle correlations. In order to take non-perturbatively at least an infinite subset of diagrams of {Q m ) into 

account, one could follow an idea of Phillips and Wallace |^,^,||4j, who investigated the three-dimensional reduction 

of the two-fermion Bethe-Salpeter equation with general four-dimensional (i.e. non- instantaneous) interaction kernels. 

Their method provides a generalization of a former formalism of Klein [^5|[36| using the quasi-potential approach of 

Logunov and Tavkhelidze |29 and has a close connection to standard time-ordered perturbation theory. Applying 

their idea to our three-body case, their method essentially consists in a systematical prescription to determine order- 

(2) . ... (2) 

by-order (in the coupling of Kj^ ) an instantaneous three-particle irreducible kernel Kj^f ^^^^ , where irreducibility is 

now defined with respect to the Salpeter propagator (Gqa/), such that 



(Gm) = (Com) - i (Com) -^^m inst(^M )• 



(132) 



This would allow for a formulation of the Salpeter equation (128) in a form that is the same as in the previously 
discussed case of vanishing two body interactions, i.e. 



<Pm = -i (Gi 



OM) 



K 



M inst 



^(3)" 



J\/- 



(133) 



We should note here already that t he Salpeter equation generally possesses both positive and negative mass solution due to 
the CPT-symmetry (see subsect. 5.2.3). In this respect both cases M > and M < have to be considered. 
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However, the method of has an inconsistency pointed out by the authors themselves: Obviously, eqs. (132) and 

( |133| ) are in clear contradiction to the occurrence of mixed energy components discussed above: due to the projector 

property of {Gqm) 



(2) 

eq. (132) would restrict {Gm) 



and consequently (Pm to purely positive and negative components 
only. In the next subsection we will therefore improve our reduction procedure such that this method nevertheless 
becomes applicable without revealing such inconsistencies. 



4.4.2 Reduction to a Salpeter equation in Hamiltonian form 



In this subsection we present a systematical reduction procedure, which avoids the difficulties and inconsistencies of the 
foreg oing first attempt and allows for a formulation of the Salpeter equation with a structure quite similar to that of 
sect. 4.3, where the dynamics was given by a connected instantaneous three-body interaction alone. Furthermore, this 
procedure will provide a systematic approximation of the exact reduced equation that is still manageable in practice 
and appropriate for explicit calculations. Our aim is to get a reduction of the Bcthe-Salpeter equation which even in 
the presence of unconnected two-quark kernels exhibits the same form and properties as the Salpeter equation ( pl| ) in 
the case of vanishing two-body terms. Consequently it then can likewise be formulated as an ei genvalue problem (or 
at least as a generalized eigenvalue problem) in Hamiltonian form as discussed in subsect. 4.3.2. In other words: 



The free three-quark propagation shall be given by the Salpeter propagator (Go a/)- Accordingly, we search for an 
instantaneous three-particle irreducible kernel V^f (a quasi potential) with irrcducibility defined with respect to 
the propagator (Goa/), which covers all the complexity arising from the unconnected two-particle interactions and 
adds to the genuine instantaneous three-quark kernel V^^'^^. 
Due to the projector property of the Salpeter propagator 



A{Gom) — {Gqm) A — (Gom) 



(134) 



where A and A are the Salpeter projectors defined in eq. ( |lOC| ) and (101), the reduced equation then merely involves 
the purely positive and purely negative energy components. Consequently the reduced amplitudes emerging from 
the Salpeter equation and its adjoint must be eigenstates of the S alpet er projector A and A, respectively. 



— However, as demonstrated in the previous discussion of subsect. 4.4.1 



the Salpeter amplitude itself is no longer 
an eigenstate of the Salpeter projector when two-particle interactions are taken into account. We found that in 
this case also the mixed energy components occur. Consequently, the reduced amplitude resulting from our desired 
reduced equation can not be the full Salpeter amplitude 'Pm but only its projected part := A<1>m- To summarize, 
we are looking for a reduction of the Bethe-Salpeter equation of the form 



m - (Go 



Ml 



1/(3) + v. 



M 



M 



with (Plj := A<P 



(135) 



Equivalence to the Salpeter equation (128) then requires that all interactions via the mixed components must be 



effectively taken into account in the quasi-potential V^f and moreover there must be an interconnection which 
allows to regain the full Salpeter amplitude <1>m and finally the full Bethe-Salpeter amplitude xm from the projected 
amplitude (P^j. 

Now let us become specific and show how such a kind of reduction can indeed be achieved. To this end we split the 
instantaneous three-body kernel V^^^ according to 

T/(3) =y(3)+yCT^ (136) 
with V^^"* that part of l^^^) which couples exclusively to purely positive and purely negative energy states, i.e. 

vf{p^,p^; p^,p;) :=:i(P5,p,) v^('Hp?,Pr,; p^,p;)vi(p^,p;) (137) 

and the residual part := V^^^ — vj(^\ which couples also to the mixed energy components and has the property 

A V^^^ A = 0. (138) 
Then we have for the Bethe-Salpeter equation and its adjoint: 



XM 

Xm 



-i Go 



M 



V, 



(3) 



(3) 



K 



Xm, 



~1 XAf M + ^^R +^Af <-^0A/- 



(139) 
(140) 
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Recall that in the case of vanishing two-particle kernels only the first part V\ contributes to the Salpeter equation, 
while the residual part V^"^ disappears according to property ( |13§| ), as discussed in subsect. 4.3.1. But now, in 

(3") . 

connection with retardation effects of the unconnected two-particle terms, also the residual p art Vp gives contributions 
to the reduction of the Bethe-Salpeter equation as has been shown in the pre vious subsect. 4.4.1. Keeping this in mind 
we now want to proceed in a way similar to our first attempt in subsect. 4.4.1 , where we transfered the effects of 

(2) (2) 

the retarded two-particle terms Kj^f into the Green's function Q^j . However, our discussion indicates that it is even 

more convenient to absorb together with K^j also the instantaneous kernel into a Green's function Qm, since the 

contributions of V^^^ occur exclusively in connection with k\^j . In this way we achieve that really all complications 
resulting from the unconnected two body-terms are gathered in the resolvent Qm^ which now is defined by 



^0 M 



(3) 



iK 



M 



1 + 1 



Qm = I 



and thus is the solution of the inhomogeneous integral equations 



Gqm — i Go 



M 



M 



Go A/ - i Q 



M 



(3) 



G, 



OA/- 



With Qm the Bethe-Salpeter equation and its adjoint can be rewritten as before in the equivalent form 



Xm 
Xm 



Qm V'f'^ Xm, 
'iXA/ Qm, 



(141) 

(142) 

(143) 
(144) 



which is suited for the six-dimensional reduction, because the new three-body kernel V)^ is instantaneous. The 
reduction is performed as before. Similar to eq. (|83|) we obtain first 



1/(3) 
va Xm 



Xm Va^ iPi,Pn) 



-A3} 
•^A 



(P«:P»)) 
(P«,P.)), 



(145) 
(146) 



where the Salpeter amplitudes <Pm and are the reductions of the corresponding Bethe-Salpeter amplitudes as 
defined in eqs. (|^ and (^6|). Inserting this back into the Bethe-Salpeter equations we get the prescription how to 
reconstruct the full eight-dimensional Bethe-Salpeter amplitudes from the Salpeter amplitudes: 



XM 



'I Qm ^M, 



XAf = -i ^Af VX Qm- 



However, with the definition V^'^'' = AV^^^ A, we in fact have 



Xm 



= -i Qm AV^^^ <!>{ 



XM^--^^iv^^^Ag_ 



jm 



(147) 
(148) 

(149) 
(150) 



showing that for a reconstruction of the full eight-dimensional Bethe-Salpeter amplitudes it indeed suffices to know 
only the projected components 



^m (p«>P»?) ^(PC^Pr,) ^M (PcP,,) , 



(151) 

(P«,P^) :=^^A/(P5,Pr,) ^(PC,P„) (152) 
of the Salpeter amplitudes. Performing now the integration over p^^ and on both sides of eqs. ( |149| ) and ( p_50[ ). 



the Bethe-Salpeter amplitudes xm and Xm the left hand side reduce to the Salpeter amplitudes ^m and <?Af and 
on the right hand side we obtain the reduction {Qm) of the resolvent Qm leading to the interconnection between the 

full Salpeter amplitudes <?m and <Pm and their corresponding projected parts <l>j^i and 'Pm, respectively, i.e. 



^m = -i {Qm) AV^''^ 
^m = -i^tiV'^^U{QM). 



(153) 
(154) 
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Here the mixed energy components of the full amplitudes <Pm and 'Pm reenter via the mixed energy components of 
{Qm)- To get the equations for the components <Pj^,j and finally have to perform the projection on purely 

posit ive a nd purely negative energy components via the Salpeter projectors A and A on both sides of the eqs. ( 153 ) 
and (154), respectively. We then find 



Ml A 



(155) 



(156) 



(157) 



where we introduced the symbol {Qm) A to denote the corresponding projection on {Qm), 

{Qm)a ■■= A{Qm)A, 

which cuts off the mixed energy components on both sides of {Qm)- Thus, due to the Neumann series of Qm, the 
reduced propagator {Qm)a may be represented as power series which starts in lowest order with the free Salpeter 
propagator (Com) and consists of an infinite number of reduced Feynman diagrams, which all are restricted to purely 
positive and negative energy components, as the Salpeter propagator (Gqa/) itself: 



{G- 



MIA 



(G, 



OMI 



A 



K 



Go M \ A 



(158) 



+ ^(Gqm (-i) 



(3) 



K 



GoM (-i) 



(3) 



K 



GoM ) A + 



The idea is now to classify the reducible and irreducible diagrams in this infinite reduced series in the same way as done 
in sect. H, where the quantum field theoretical six-point Green's function G was non-perturbatively constructed as a 
solution of an inhomogeneous eight-dimensional integral equation. But now this classification is done on the reduced 
level where irreducibility is understood with respect to the 'free' Salpeter propagator (Gom)- This means that we are 
looking for an irreducible kernel V^f such that {Qm)a is the solution of the following six-dimensional integral equation: 



{G 



MA 



(Go 



Ml 



i (Com) Vm {Gm)a = (Go a/) - i {Qm)a V^f {Gqm). 



(159) 



Note that in contrast to our first attempt, where this ansatz due to the restrictive action of the Salpeter propagator 
{Gqm) on purely positive and negative energy components led to inconsistencies, here the ansatz becomes possible 
now, because {Qm)a itself and thus all terms of the series (15S) have by construction the same restriction as (Gom) 



to these components only. Formally the determination of corresponds to the inversion of {Qm)a, which due to 
the projector properties is restricted to the subspace of positive and negative energy components. In particular, this 
requires the inversion of the Salpeter propagator (Gojv/) ii^ this subspace. For this purpose we introduce the operator 
hoM by 



-1 7 



'7 



.7° [A/ I - Ho (PC, p,,)] 



X (271)3^(3) (p^^_p,) (2^)3^(3)(p^_p;), 



(160) 



with Ho the free three-fermion Haniiltonian defined in eq. (110) such that the 'inverse' of (Go a/) in this subspace is 
given by^ _ 

(GoAf) hoM = A, hoM (GoA/> = A (161) 
and the 'inversion' of {Qm)a can now be expressed by 



{Qm)a [ hoM + i V^f ] ^ A, [ hoM + i ] {5m) a = ^ 



(162) 



A unique definition of the effective, irreducible kernel V^f then requires its restriction to positive and negative com- 
ponents according to 



AV, 



cff 



V, 



off 
M • 



(163) 



Finally, having found this quasi potential VJ 



off 



M 

, we can use eqs. (162) in order to transform the Salp eter equation 



(155) and its adjoint (156) into the desired form as indicated in the beginning of this subsection in eq. (135) 



M - {Go 



Ml 



y(3) -I- v; 



M, 



(164) 



"^M = -1 ^M 



1/(3) + v; 



M 



{Go 



Ml 



(165) 



® Note that according to our concise operator notation A and A here have the meaning of an integral operator i.e.: 

^(Pe,P.; P?,p;) :=^(pe,p^) (2^)^ <s(^)(p« - p?) (27r)3 5(3)(p, - p^) 
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The form and therefore the properties of this reduced bo und state equation are indeed exactly the same as in the 
case of vanishing two-particle kernels discussed in the sect. 4.3. The only extension is the effective quasi potential 
which occurs in addition to the genuine instantaneous three-body kernel V^^"^ and absorbs all the complexities entering 
via the retardation effects from the unconnected two-quark interactions. Note, however, that V^f , in contrast to V^-^\ 
in general is energy-, i.e. M-dependent as indicated by the subscript M. 



To obtain a Hamiltonian formulation of the Salpeter equation, we multiply eq. (164) by [ i 7° (g) 7° (g) 7° /iqm ] 
and thus end up with a the generalized^ eigenvalue problem which determines the bound state mass M and the 
corresponding amplitude <?^: 

(166) 



Mi 



where the Salpeter Hamiltonian TLm now explicitly reads 

d'p^ d^p' 



+ ^(Pc,P„) 7°'»7°'^7'' 



(27r)3 (27r)3 
d'p^ d^p; 



(27r)3 (27r)3 



+ ^(P?,P„) 7°«)7°®7° 
The next step is to determine V^f . According to eq. (157) we have 

{Qm)a - A{gM)A, 



^^^'Hp«,p.; p^,p;X(p^,p;) 



(167) 



(168) 



where on the left {Qm)a is given by the integral equation (159) which defines V^f and on the right we insert Qm as 

given by the integral equation ( |142| ) with kernel + if . This equation then has to be solved for V^f . As shown 
in detail in appendix [A| with the restriction (163) the effective instantaneous kernel V^f can be uniquely determined 
order-by-order as an infinite power series expansion 



off 
M 



(fc) 



fe=i 



of irreducible interaction terms V^f'^^'^ in powers k of the kernel V-^' -\- K\J . The explicit expressions in arbitrary 
order k are then constructed according to the following prescription 



(169) 



(3) , 1^(2) 



v; 



=ff(i) 



M 



hoM A ( Gqm Gqm ) A /iqmj 



(170) 
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i A ( Gqm (-i) 



(3) 



K 



G, 



OM 



(-i) 



y(3) ,7^(2) 



k til 



Gqm ) A h 



QM 



E 

^■-2 ki, k2, . . . ,kr < k 
ki + k2 + ■ . ■ + kr = k 



"1 



{Go 



Ml 



"1 ^M 



{Go 



Ml 



1 ^M 



(171) 



Notice the emerging structure of these equations: The reduced Feynman diagram of fcth order {i.e. the first term on 

the right hand side of eq. ( |171| )) consists on the one hand of the irreducible part Vm'''^^ of order k, which we are in fact 
interested in, and on the other hand it contains an order- A; reducible part, built from all possible iterations of reducible 

diagrams V^^'^'^^ of lower order ki < k with ^ ■ ki — k, as given by the second term in ( |l7l[ ), which obviously has to 

be subtracted to get the desired Vj 



=ff(fe) 



M 



'generalized' means that now the Salpeter Hamiltonian itself depends on the eigenvalue M due to the -dependence 
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4.4.3 The normalization condition for the reduced amplitudes 



The solutions Xp of the Bethe-Salpeter equation (12C) have to satisfy the normalization condition, which may be 
formulated in the explicitly covariant form (|69[). The covariant framework ensures that the normalization in the rest- 
frame implies the correct normalization of the Bethe-Salpeter amplitudes Xp in ^-^y frame. In this section we will 
determine the correspo ndin g normalization condition for the projected Salpeter amplitudes i-e. the solutions of 
the Salpeter equation ( 164 ) in the rest- frame. For this purpose we start with the Bethe-Salpeter norm which in the 
center-of-mass frame reads 



1 Xm 



XM 



2M, 



(172) 



where the pseudo-Hamiltonian Hm is defined by 



Hm — G, 



om 



(173) 



To get the analogous condition for the reduced amplitudes, we have to express the eight-dimensional Bethe-Salpeter 

amplitude xm and its adjoint Xm by the corresponding reduced six-dimensional amplitudes <Pj^j and ^jv/, respectively. 
We do this by using the relations (149) and (150), i.e. the prescription how to reconstruct the Bethe-Salpeter amplitudes 
from the Salpeter amplitudes: 



XM 



-i Qm AV'^^^ $i 



Mj 



(174) 



Xm 



M- 



(175) 



Recall that in eq. (141) the Green's function Qm was defined as the resolvent of the pseudo-Hamiltonian 



rrR 



I.e. Hit Gm = Gm Hft = I. 



Accordingly, with the decomposition V^^'^ = V^^^* -I- , we write 



.(3) 



Hm — Hfi + i V"j^^ , 



(176) 



(177) 



where the projected part V\ of the instantaneous kern el V^ '^^ has no explicit M-dependence as V^^'^ itself. Therefore, 
it gives no contribution to the normalization condition (172) which thus becomes 

d 



i^i V^^^A G 



M 



rrR 



Gm a 



M 



2M. 



(178) 



Using the resolvent equation (176) for the Green's function Gm the derivative of Hfj can be rewritten as a derivative 
of Gm- 



d 
DM 



H 



m 



Gm — 



d 
"dM 



^m- 



(179) 



Substitution into (17S) then yields 



i KV^'^A 



dM 



JM 



A V^(3) = 2M, 



which is convenient for the further reduction, since the 'vertex' amplitudes V'^'^^A and A F^^^ <?;^j are six- 
dimensional and the only remaining full eight-dimensional quantity is the derivative of the Green's function Gm which 
can be reduced by the internal integrations over and according to eq. (p3[): 



(180) 



(181) 



The next step is to get rid of the three-body kernel V^^'' by absorbing it into the Salpeter amplitudes 'P^ and <P]^.j by 
means of the Salpeter equations (155) and (|156|). We start with 



(i^Sm)a = A [a{Gm)a] = A (g 



\dM 



dM 



m)a ' 



(182) 
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where we used that the partial derivative d/dM can be c omm uted with the integrations over p^'^ and p^*' and also with 
the Salpeter projector A. With the resolvent equation (162) for the reduced propagator {Gm)a, we can then rewrite 
the derivative of {Gm)a as 



d 



{5m) A [hoM+iV^f]{g_ 



(183) 



which finally enables us to apply the Salpeter equation (155) and its adjoint ( |l56| ) in order to absorb the three-body 
kernel V^'^'^ into the Salpeter amplitudes: 



i < V-''^ {Qm)a irr [ + i ] {Qm)a 1^^'' < = 2M 



dM 



dM 

The explicit expression for the derivative of the operator /iqmj owing to its definition (|16C|), is given by 



d 
'dM 



'-0M 



(PCP,; P^,p;) = -i7°»7''®7° (2^)3 5(3) (p^_p,) (2^)3 5(3) (p^_p.^)^ 



(184) 



(185) 



and one readily shows that the general relation (111) between the rest-frame Salpeter amplitude and its corre- 



sponding adjoint holds likewise for the projected amplitudes = A^m and = ^mA, i.e. 



-tA 



<?M(P^Pr,) =-<^M {Pi,Pn)T®T^T^ 
so that we finally end up with the following form of the normalization condition for the reduced amplitudes <l>^j 



(186) 



2M. 



(187) 



Here (-j-) denotes the positive definite scalar product (117), which is induced by the /^^-normalization condition of 
the Salpeter amplitude in the case of vanishing two-particle kernels. In comparison to the case where the dynamics is 
determined by an instantaneous three body kernel alone we thus find that (owing to its explicit energy dependence) 
the effective kernel V^f in general leads to an additional contribution to the norm. 



4.4.4 Lowest order contributions to the effective kernel 



In equations ( |l69| - 171) we displayed the general order-by-order prescription to construct V^f . In practice, we have to 



approximate the effective kernel V^,f , which itself consists of an infinite number of terms. A systematical approximation 
is now given by truncating the series ( |169| ) at some finite order fc < oo, i.e. 

< V,f^'' + V^,f^ +... + Vif '\ (188) 

thus yielding an approximation of the Salpeter amplitude c:i by the solution of 

= -i(GoM) {v^'' + E V^f'^ <\ (189) 

Note that such a finite order approximation of V^f means for the original reduced propagator {Qm) A (and thus also 
for the Salpeter equation) an approximation beyond perturbation theory, due to the infinite iteration of V^f in {Qm) A- 
It is worthwhile to mention that this subsequent approximation of the Salpeter equation (within the CMS frame) still 
preserves formal covariance. Our reduction procedure and thus the construction of V^f can be covariantly formulated 
in any arbitrary reference frame according to the covariant replacements — > pn and p p± as mentioned previously. 



Consequently, also the truncation (188) of t he ef fective kernel V^f can in fact be performed frame- independently. An- 



other aspect concerning the approximation (188) requires attention. For the calculation of transition matrix elements 
we need the full Bethe-Salpeter amplitude xm which (if and (l>^,f are known exactly) can be reconstructed by the 

prescription ( |l49| ) via the Green's function Qm- To be consistent we need an approximation Q^j^"^ that corresponds to 



U. Loring et ai: Relativistic quark models of baryons with instantaneous forces 



33 



the approximation (188) of the effective kernel In other words, we require the corresponding order k approximation 



Xm of Bethe-Salpeter amphtude x m such that its reduction yields the order k approximation <P^j ^ of the Salpeter 
amplitude. As shown in ref. |3l[] a consistent prescription for an approximated reconstruction of the Bethe-Salpeter 
amplitude can indeed be found. 



With regard to explicit calculations let us now become specific and compute the explicit expressions for the 
contributions to V^f up to second order. 



The Born term V^f 



(1) 



Concerning the Born term V^f''^^ we can refer to a former result of subsect. 4.4.1 given in eq. (131). According to eq. 
(11701) the projectors A and A select the purely positive and negative energy components, i.e. the terms (a) and (b) of 



eq. ( |13l| ), and cut off the mixed contributions, which correspond to the terms (c) to (f). We then find the following 
rather simple result, which is a sum of three unconnected two-fermion potentials for each quark pair (see also fig. 
for the diagrammatic representation): 



vTf^'^(p«,P,; Pop;) = 



hoM A { Gqm Gom) a hoM 



7° (g) 7° (g) 7° X 



A^ (g) A. 



2 A^ 



7"® 7° F(2)(p^,p^)l ® I (2^)3 J(3)(p^-p;) A+'^A+'®A+' 



-A^ (gylj 



7°®7° F(2)(p^,p^)j ® I (27r)3 J(3)(p^_p;^) A-^' ® A^' ® A^' 



cyclic, perm, of (12) 3 corresponding to the 
interacting quark pairs (23) and (31) 



(190) 



where ylf = /If (pi) and A^ = /if (pj)- Note that this Born term V^f^^' in fact is M-independent 



ff(i) 



H 

3+- 



I 



I 

3'- 



— 2 



— 3 



Fig. 12. Time ordered graphs for the Born term {Gom) Vjv*^^' (Gom) due to eq. 
(shown for the quark pair (12) only) is represented by the vertical dashed line. 



(19C). The instantaneous two-body kernel 



The second order term Vl 



=ff(2) 



M 



Already in second order the expressions become much more complex: Using the recipe given by eqs. (169 ~ 171) we 



obtain for the second order term V^f^^'': 



v; 



=ff(2) 



M 



-i I h^M A (Gom K^m Gom) (I - A) V<^^^A + A V^^\l - A) (Gqm K^S Gom) A hoM 



*0M 



A (^GoM -^M GoM Gqm) a 
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^ ( Com Com ) ^ A ( Gqm -^m 



A 



ho 



AI 



(191) 



Analyzing eq. (|191[) in more detail, V^/*'^^ essentially consists of three structurally different contributions, 



v; 



=ff(2) 



M 



(2) 
M • 



(192) 



The term w'^^j is given by the first term on the right hand side of eq. ( |l9l| 



(2) 



M 



i /loA/^ (Gom Com) (I - A) T/^^) yi 



i ^ (I - yl) (GoAf ifi/ GoA/) /I /io 



(193) 

and again we can go back to 
subsect. 4.4.1 and use the result (131) to comput e th e explicit expression of Wj^/ . However, in contrast to the Born 



It is of first order in the instantaneous two-particle kernel (as the Born term VJ 



cff(l)^ 
M ) 



term, now the terms (c) to (f) of the expression (131) that couple also to the mixed energy components, enter only. 



They are attached in a symmetrical way from the right and left hand side to the residual part vS^^ of V^^\ such that 



(2,') 

the mixed energy components appear internally in Wl/ . Notice that the terms of W^^ are suppressed with respect 



(2) 



to the corresponding reducible terms vj(^\GoM)yM 

off(l) 



' M 
(1) 



T/off(l); 



{Gom)Va^\ built by iteration of V^^' and the first 

order Born term VJ/*^ '' by means of the Salpeter equation. This is apparent from the different non-singular internal 

ener gy d enominators in Wj^,/ in comparison to the singular expression (Go a/) that emerge in the reducible t erms (see 
eq. (131) and the subsequent discussion of the corresponding terms (a), (b) ^ (c), (d), (e), (f) in subsect. 4.4.1 ). A 
diagrammatic representation of the corresponding time-ordered Feynman graphs is shown in fig. n3. 



(3) 



2^ 



+ 2- 



t3' 3' 




Fig. 13. Time ordered graphs of the contribution to the second order term Vm'^'^^ given in eq. (193). The instantaneous 

three-body kernel is represented by the vertical dashed line that connects three quark lines (indicated by the dots). The 
instantaneous two-particle kernel is shown for the pair (12) only. 



Furthermore, the second irreducible part on the right hand side of eq. (191), which is of second order in the two-body 

kernel, may be split into two terms of different structure. According to the decomposition K^^ — K^^ + Kf^ + Kf^, 
with 



Kii{Pi,P,v PlPn)^V^'Hp^,Pi)^Sl '(iM-p,) (2^)4<5(4)(p,,_p;j 



(194) 



and iffl , Kf} the corresponding cyclic permutations of Kl^ , firstly we find an unconnected part C/lJ-* that consists of 
a sum of irreducible two-body loops in each quark pair, i.e. 



U 



A h, 



OM 



A (Gqm Go a/ -^aI Go a/) A 

- A (Gom Klj Gom) a hoM A (Gqm ^a/ Gqm) A 



(195) 



h. 



OM 



corresponding terms with interacting quark pairs (23) and (31). 
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Remember that the second term in eq. (195) just subtracts the reducibl e pa rt of A {Gom -^1/ Gqm K\j Gqm) A that 
is built up by a two-fold iteration of the corresponding Born graphs (19C), so that we are left with an irreducible 
double Z-loop graph of the corresponding time-ordered Feynman diagram as shown in fig. Q4[ 



Fig. 14. Time ordered graphs of the unconnected irreducible two-particle kernel 11^^ as defined in eq. (195). Here only the 
term (12)3 is shown. 



(2) 

Secondly, we find a connected part C\/ that contains the sum of all possible irreducible quark-exchange diagrams. 



This term is given by 



A (Gom Gom K'^j Gqm) ^ 
— A (Gom Klj Gqm) ^ ^om ^ (Gom Gom) ^ 



(196) 



+ 



corresponding terms with other quark pairings in the incoming and 
outgoing channels: (23, 12), (12, 31), (31, 12), (31, 23) and (23, 31) 



In fig. 15 the different time-ordered graphs contributing to the irreducible second order quark-exchange interaction 
are shown diagrammatically. 
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3 + 




(2) u ^ 

Fig. 15. Time ordered graphs of the connected irreducible quark-exchange interaction C]^/ as defined in eq. (19f). As an 
example only term (12, 23) is shown. 



The explicit calculation of the second order terms yl(GoM ^om -^^m ^om)A, needed for the determination 

(2) (2) 

of Ulf' and C}^' , is lengthy but straightforward and can be performed by making again elaborate use of the residue 
theorem. Owing to the increasing number of quark-lines and the increasing number of contributing poles in the relative 
energy variables andp°, the structure and coordinate dependence of these explicit expressions is rather complicated 



in comparison to the rather simple structure of the Born term (19C). Moreover [/{^^ and C)^/ exhibit an explicit 



M-dependence. We restrict our explicit calculations to the leading Born term (Born approximation) 



5 Bound-states in Born approximation of the quasi-potential 

The discussion of the lowest order contributions to the effective quasi-potential V^f in the last section clearly showed 
that with increasing order of the contributions to V^f the explicit expressions rapidly become more complicated. While 



36 



U. Loring et ai: Relativistic quark models of baryons with instantaneous forces 



the dominant leading order Born term V^f is still rather simple in structure, already the second order contribution 

^M^^'' contains a lot of different irreducible terms whose structure is quite complex and thus impedes an efficient 
numerical treatment. Therefore, expecting these contributions to be small in comparison to the leading Born term, we 
consider the Born approximation 

:^ V^f'^ (197) 

only. For the sake of completeness let us summarize the corresponding expressions for the Salpeter equation and 
the normalization condition in this approximation. These equations shall constitute the basis of our quark model for 
baryons. 



5.1 Salpeter equation and normalization condition 



Approximating the series ( 169 ) by the leading Born term ( 190 ), the approximated Salpeter equation ( 166 ) can still be 
formulated as an ordinary eigenvalue problem. 



M 



M <P 



(198) 



since, due to eq. ( |l90|) , the Born term V^f^^"^ in fact is Af -independent. The M-independent Salpeter Hamiltonian Ti. 
then reads explicitly: 



[n<Pii] (PC, p.,) = no{Pi,Pn) K{Pi,Pv) 

+ [A+{pi) ® A+{p2) ® A+{p3) + /Ir(Pi) ® A^{p2) ® (P3)] 



(199) 



+ [/If(pi) <» A+{p2) ® A+ips) - A^ (pi) ® A^ (pa) ® A^ (pg)] 

X 7%7°®I J ^V^^Hp^,p'^)^l<Pij{p'^,p,,) 
+ corresponding terms with interacting quark pairs (23) and (31). 



Note the striking structural difference between the connected three-body part and the unconnected two-body part: 
The two body term shows a relative sign between the positive and negative energy projectors and occurrence of the 
identity (instead of 7°) in the Dirac space of the spectator quark. 



To be consistent, the same approximation of V^f must also be us ed in the normalization condition (187). In Born 
approximation, the second term in the normalization condition (187) vanishes, owing to the explicit Af- independence 

of the Born term V^f ^ and we arrive at 



V 7 



dM 



M 







2M. 



(200) 



Consequently, the solutions <I>f.j of the Salpeter equation ( 199 ) in Born approximation of V^f have to fulfill the same 
>C^-normalization condition (11(:) as in the case where the dynamics was determined by the instantaneous three-body 
kernel alone (see subsect. 4.3.3), i.e. 



(201) 



The Salpeter equation for <Ppj and the corresponding adjoint equation for ^j^j must be consistent with the relation 



( |186| ) between and 'P'^j. This leads to the following condition for the instantaneous interaction kernels V'^'^^ and 



7" «) 7° «) 7° 



it 



v^'\p'^,p',, Pe,p,) 7"®7°®7" = ^^'Hp?,p„; p?,?;,) 



7° 7° 



V^'\p',,Pi) 



y' = v^'Hp,,p',h 



(202) 
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which impUes that the Salpeter Hamihonian (199) in Born approximation of the effective kernel is hermitean with 
respect to the scalar product (-I-), i.e. 

{$i\n $2) = {n <Pi\^2) V <Pi,<l>2 with ^<?i^2=<?i,2 (203) 
As in the case of vanishing two-quark kernels this again guarantees that 
~ the eigenvalues (bound-state masses) M oiTi, are real, i.e. M* = M; 

— the Salpeter amplitudes and corresponding to different eigenvalues Mi ^ M2 are mutually orthogonal: 
(^11*1) =0. 



5.2 Symmetries of the Salpeter equation 

So far we discussed the constraints ( ^02|) of the instantaneous two- and three-quark interaction kernels y^^-* and ^^^^ 

that followed from the interconnection of the amplitude <?^j and its adjoint <?^j and guarantee the hermiticity of the 
Salpeter Hamiltonian Ti, with respect to the positive definite scalar product (-I-). We are led to further conditions on 
the kernels if we regard the symmetries which the strong interaction of the quarks has to respect. Since the underlying 
theory, quantum chromodynamics (QCD), is invariant under parity transformations {V), time-reversal (T) and charge 
conjugation (C), these symmetry properties must be incorporated in the Salpeter equation. This means specifically: If 
is a solution of the Salpeter equation, the same must also hold for V with V g {7^,T,C} the representation 
of the corresponding transformation on the (projected) Salpeter amplitudes = A(I>m- Below we shall investigate 
the corresponding constraints on the interaction kernels V^^^^ and T^^^-* that follow from these invariance conditions. 
Instead oiV, T and C we alternatively consider V, T and CVT. 



5.2.1 Parity invariance 



The representation of the parity transformation 'P{x'^j x) := {x^, — x) on the full momentum space Salpeter amplitudes 
<?A/ is given by 

[T'<I>m](pc,P^) =7°®7°®7"'^A/(-Pc,-p„). (204) 

Owing to the intertwining relation yl^(pi)7'^ = 7*^7!^ (—pi) the different energy components of <?m represent invariant 
subspaces under the parity transformation V, such that V decomposes into irreducible representations on these different 
subspaces. In particular, the Salpeter projector A commutes[^ with V, i.e [P, yl] = such that 

[^<?m] (PCP^) =7°®7°®7°'Z'm(-Pc--P^) with AVK = VK (205) 

is the rep resen tation of V on the projected Salpeter amplitudes <P^j = A(1>m, which actually appear in the Salpeter 
equation (198). Parity invariance implies that with also V^P^j is a solution of the Salpeter equation, i.e. the 
Salpeter Hamiltonian has to commute with the representation V of the parity transformation, i.e. [V, Ti.] = 0. With 
[P, yl] = and the invariance of the free Hamiltonian TCq under V, i.e [V, Ho] = 0, one readily deduces the following 
conditions for the three- and two-quark interaction kernels: 

7°®7"®7" W3)(-p5,-p,; -p^,-p;)7"®7°®7" = ^^'^(PC,P.; P'^P'^)^ 

7°®7° ^^'H-P6-Pc) 7°®7° = V^^Hp^,p'^). (206) 

As usual, [P, H] = also implies that the solutions of the Salpeter equation simultaneously are eigenstates of V, 
i.e. 

V-Pi^^ = 7T^ij^^, (207) 

with definite parity n — ±1. 



The brackets [■, ■] denote the commutator [A, B] := AB — BA 
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5.2.2 Time-reversal invariance 

The representation of the time-reversal transformation T(x'^,x) :— {~x'^,x) on the full momentum space Salpeter 
amplitudes reads 

[T^m] (P6 P„) = -7^7' ® 7S' ® 7S' <?m(-Pc, -P^). (208) 

Again the different energy components of $m define invariant subspaces under the time-reversal transformation ac- 
cording to the intertwining relation A^{pi) 7^7^^ — 7^7"^ Af (— Pi). In particular, we find that the Salpeter projector 
is time-reversal invariant, i.e. [T, yl] = 0, such that we have a representation of the time-reversal transformation on 
the subspace of purely positive and negative components, i.e. for the projected amplitudes 'Pf.j — A<1>m holds: 

[r<] (P5, p,) = ® ^ ^ly <*(-p^, -p,) with A r<?f, = T<Pij. (209) 

To respect time-reversal invariance of the strong interaction, we must impose that [T, Ti.] = 0. Using the invariance 
property [T, vl] = of vl and the time-reversal invariance of the free Hamiltonian Ho, i.e. [T, Ho] = we end up 
with the conditions 

- 7'7^^7'7^^7'7^ V^^'^X-p^, -p^; -p^, -p^) j^-f^®-f^-f^®j^j^ ^ V^^) (p^^ p^. p^, p^), 

7V ® 7^7' V<^^^\-p^, -p^) 7V ® 7^7' - V<^^\p^,p'^). (210) 



5.2.3 CT'T-symmetry - Interpretation of negative bound-state masses 



The Salpeter Hamiltonian H being hermitean with respect to the positive definite scalar product ( |117[ ) guarantees that 
the eigenvalues M, i.e. the bound-state masses, are real, as one imposes for physically acceptable solutions. However, 
H is not positive definite, since even the free Hamiltonian Ho is not positive. Accordingly, H in general possesses both 
positive and negative eigenvalues and the spectrum might be even unbound from below. These negative eigen value s, at 
first face seem physically unacceptable and the corresponding amplitudes also contradict the normaUzation (201) via 
the positive definite £^-norm. Nevertheless, these negative energy solutions with M < can be interpreted physically 
meaningful. In fact, since our covariant Salpeter approach is based on relativistic quantum field theory, it should reveal 
a particle-antiparticle symmetry as a characteristic feature due to CPT-invariance. Accordingly, we demand that the 
instantaneous two- and three-quark interaction kernels commute with the Dirac-space operator ^^^^i j'^ "f^ , i-S- 



7V' 



'7V 



V^^\pi,p„: p^,p;) 
7V® 7V, V^'\p^,p'^) 



= 
= 



(211) 
(212) 



in order to ensure that the three-quark Salpeter equation in fact respects the CPT-symmetry of the strong interaction. 
In this manner the negative bou nd-st ate m asses get a well defined physical interpretation as we will see in the follo wing 
discussion. The conditions (211) and (|212|) on V'^^'' and V^^'' imply that the Salpeter Hamiltonian H given in eq. (199) 



anticommutes with {^^^x 7*^7^ 



{7°7^ 



.7V, n} = o 



as can easily be shown with the anticommutatorp] and intertwining relations 



{7V, 7°} = 0, {7V, H^iPi)} - and 7V A*(Pi) = ^?(Pi) 7V 



(213) 



(214) 



Moreover, it follows from ( ^14[ ) that 7" 
(204) anticommutes with 0i^i7"7^: 



7" and hence also the representation V of the parity transformation 



, 7O75 : 



7V, V} = 0. 



(215) 



Now let 'I'^j^j TT be a solution of the Salpeter equation with negative mass —M < and parity tt which obeys: 



H $ 



A 

-M, 



A 

-Af. 



(216) 



We consider the transformation <Pi 



M, TT 



^^j^j „ of the amplitude <l>^]^f ^, given by 



^-M,^(P«,Pr,) 



7V' 



7V <?^!m,.(p«,p.)- 



(217) 



The brackets {•, •} denote the anticommutator {A, B} :— AB + BA 
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Then, due to eq. (213), also this CPT-transformed amplitude rr is a solution of the Salpeter equation, but now 



with the positive bound-state mass +M > 0. At the same time eq. (215) implies that parity — tt opposite 



to 'I'^j^f ^. Thus, we have 



M, 



= +M $i 



M. 



M, 



M, 



(218) 



Consequently, the eigenvalues come in pairs with opposite sign, but with eigenfunctions (Salpeter amplitudes) having 
opposite parity. This symmetry indeed allows the interpretation of the negative energy solutions for a given set 



of quantum numbers as antibaryon states, which after the transformation <P"^ 
solutions of opposite parity but otherwise with the same quantum numbers: 



M, TT 



'^-M TT yield positive energy 



i=l 



M, -TT- 



(219) 



This is a new interesting feature of our Salpeter equation-based baryon model in contrast to nonrelativistic (or 
relativized) quark potential models, which are usually based on the ordinary Schrodinger equation: Solving the Salpeter 
equation for fixed spin J yields at the same time both the positive and the negative parity bound state spectrum of 
the baryons, see fig. |l^ for a diagrammatical illustration of this feature. Furthermore the positive and negative parity 
states are coupled in this way and are not independent as in the ordinary nonrelativistic potential models. 




Fig. 16. Interpretation of the negative energy solutions due to the CPT-symmetry of the Salpeter equation represented by eq. 
(21£). See text for explanation. 



Notice that, owing to the intertwining relation 7^7^ yl^(pi) = Af{pi) 7^7^, the roles of the positive and negative 
energy components are interchanged by the CPT-transformation: 

3 3 
^M+;. = (8)7V <?:m:. and <P]^Z^ = (g) 7 V <^;|/t.- (220) 

4=1 1=1 

Consequently, only both subspaces of purely positive and negative energy components together (but not separately) 
define an invariant subspace under the CPT-transformation. Thus really both subspaces are necessary to get an 
(irreducible) representation of CVT. In particular the so-called reduced Salpeter equation, in which the negative 
components are a priori neglected (Tamm-Dancoff approximation), violates in general the CT'T-symmetry. 
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5.3 Eigenstates of the Salpeter projector 

Let us come back to the specific projector structure of the Salpeter equation and discuss in some more detail the 

corresponding induced structure of the solutions, i.e. of the (projected) Salpeter amplitudes (p^^ = ^^m, which ob- 
viously are eigenstates of the Salpeter projector. These projection properties of the solutions reduce the number of 
independent functions necessary to describe the baryon state. Due to the form of the Salpeter projector A the solutions 
spht into the two orthogonal purely positive and purely negative energy components and thus the (in Dirac space) 
64-component function in fact reduces to an effectively 16-component function only. To perform this reduction of the 
Salpeter amplitude to a 16-componcnt function we have to determine the general form of the eigenstates of the Salpeter 
projector A. To this end we will consider first the positive and negative energy solutions of the free Dirac equation 
for a spin-1/2 particle. These four-component Dirac spinors, which are the eigenstates of the energy projectors A^, 
can be constructed in the usual way by the embedding map of two-component Pauli spinors. This scheme can then be 
generalized to the (projected) three-fermion Salpeter amplitudes which are eigenstates of the Salpeter projector 
A and accordingly are formed by a three-fermion embedding map. 

For the following considerations it is convenient to adjust our notation to the symmetry properties of the Salpeter 

amplitude under permutations of the quarks, especially in the case of different quark masses. So far we used a simplified 
notation suppressing the flavor dependencies of the single-quark operators and Af and we assigned to each quark 
i an individual quark mass rrij. With the replacements of Hi and Af given by 

H{pi) -.= ^2 Hmf{Pi)<S>Vf, with Hmfip) ■■=1° {fP + mf) (221) 
/ 

^^(Pi):=E<.(PO^^f- with <(P):=^^^^^^^^^^^ (222) 

and (Jmf (p) := y^lp^ + mj, (223) 

the correct assignment of quark masses according to their flavor / = u,d,s is realized by the flavor projectors 
:= I/) (/I such that the free Hamiltonian Ho and the Salpeter projector 

A{p^,p„) = ^+++(p^,p„)+vl— (p^,p^), (224) 

^^^^(PcPr,) :=yl±(pi)(8yl±(p2)Ovl±(p3) (225) 
become permutationally invariant operators. 



5.3.1 Dirac spinors as embedded Pauli spinors 

Let us first consider the Dirac spinors 'ipm '■ -R-m ~^ ^'^^ ^ single spin-1/2 particle with mass m, i.e. the positive 
and negative energy solutions of the free Dirac equation on the positive/negative mass shell RJ^ := {p € ]R^ : {p,p) = 
rn?, p = {±(jJm{p)-iP)}- These are eigenstates of the positive and negative energy projectors A^: 

<(P) V'+(P) = V'+(P) with p=(+a;„(p),p)eR+ (226) 

^m(p) V-m® = V'm(p) with p = (-a;„(p), p) € R" . (227) 

As usual, the positive and negative energy solutions • R-m ^ Dirac equation may be written in the Weyl 

representation as 



^+{P) = T+{p) ^Up) with T+(p) := \ ( ^^^) (228) 

V2a;„(p) V V'^iP) J 

i'^niP) = T-{p) ^-(p) with T-(p) := ^=i== (^/M) (229) 

V2w„(p) W^i'^P)/ 

where (p^ : R^ C'^ are two-component Pauli spinors, V is the parity transformation, i.e. V p = (<^m(p)j— p) for 
P = {P°,P) = (wr„(p),p), and 

(j{p) 



a(j>):=a,i^ => V^=^^^^^, (230) 

^y2{u>m{p) + m) 
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with ai the Pauh matrices and do = cr° = Icps . We wrote these relations already in a form that defines the so-called 
embedding operations T'^(p) : t—i- C*. They map arbitrary two-component Pauli spinors ip^ into four-component 
orthogonal cigcnstates i/"™ of the energy projectors A^. This is also apparent from the properties yl^(p) T'^(p) = 
T^(p) and yl^(p) T'^(p) = of these embedding maps. On the other hand the mappings tp^ ^ (p^ are also unique, 
since they satisfy 

[T±(p)]^ r±(p) = I^. and thus ^± = r±^ ^± (231) 

and, in particular, they are isometric operations: tpm^'^'m — 'fni'fm- Finally, we define single-quark embedding oper- 
ations 

r±(p):=^ r±/p)0 7'f (232) 
/ 

which account for the correct m ass assignment for each flavor / and accordingly map to eigenstates of the energy 
projectors ^*(p) defined in eq. (222): 

yl±(p) T±(p) = r±(p) and /lT(p) T±(p) = (233) 



5.3.2 Embedding map for Salpeter amplitudes 

Now we use this result for the construction of the solutions 

<(P?,P.) ='?i/+(P?,Pr,) +<Z'Ar"(p«,pO (234) 
of the Salpeter equation, whose components 

'Z'm**(p?,P,) = ^***(PC,P,)<Z'A/(PC,P^) (235) 

are eigenstates of the positive and negative energy projectors A^^^ . According to the above di scuss ion we now define 
the three-quark embedding maps by the tensor products of single quark embedding operators (232) 



-i±±± 



(pe,p,,) :=r±(pi)®T±(p2)®r±(p3). (236) 



Then the positive and negative energy contributions ^^/^^ to can be uniquely written as 

<?M**(PC,P,) = T±±±(p^,p,)v^±(pe,p,) (237) 

in terms of the embedded three-particle amplitudes t/?^ which involve triple tensor products of Pauli spinors only. 
Finally, we have for the Salpeter amplitude the unique orthogonal decomposition 

^mCp^P')) =7'+++(PcP'7) '^MlPc.Pr,) + T—{^^,^„) v7-,(p5,p^) (238) 

and thus the determination of the (in Dirac space originally 4 ® 4 ® 4 =) 64-component solution <Pjii of the Salpeter 
equation reduces via the embedding map to finding the two (only 2 2 (g) 2 =) 8-component amplitudes <p^j- Due to 
the isometry of the embedding maps which follows from eq. (231), the normalization condition for the Salpeter 

amplitudes can be expressed in terms of these Pauli-spinors according to 

«|0 = {if+M,) + {pl.,\^l,) = 2M, (239) 

where ('/'mI'/'m) denotes the usual nonrelativistic (positive definite) norm. The next step is to investigate the 
structure of these three-quark Pauli amplitudes for a given set of quantum numbers specifying a baryon, in order to 
find a proper basis. 



5.4 General decomposition of the Salpeter amplitudes 

The three-quark Salpeter amplitude of a baryon is characterized by a set of quantum numbers that are conserved 
under the strong interaction. We consider in this work light baryons, which are built up by quarks with fiavors up (u), 
down (d) and strange (s). The fiavor-S'[/(3) symmetry is explicitly broken to SU{2)®U{1) by the different constituent 
quark masses m„ = ma < rus- With ruu = rud = rUn only the SU(2) isospin symmetry shall be assumed to be exact. 
Due to parity invariance, rotational invariance and this (broken) flavor invariance a baryon is then characterized by 
the parity tt, the total spin J with 3-component Mj, isospin T with 3-component Mt and strangeness S* . Moreover, 
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according to Pauli's principle, the (projected) Salpeter amplitude, together with its positive and negative energy com- 
ponents, must be totally antisymmetric under permutations a Cz S3. 

Consider the three-quark Salpeter amplitude = <l>^j j^MjTM s* describing a baryon with the quantum num- 
bers listed above. In order to determine the structure of its embetTded Pauli spinors, we have to investigate, how 
the corresponding transformation properties of the Salpeter amplitude transfer to the Pauli spinors ipj^ via the 
embedding maps and T : 

— The representation V of the parity transformation is given for the Salpeter amplitudes by 

(p«,p,)=7°®7°®7°'Z'M(-P«,-Pr,). (240) 



With our special choice (228) and ( |229| ) for the embedding operations T^, we find the following simple intertwining 
relations 

V 7^±±± ^ j,±±± ^ (241) 

On the right hand side, iP' is the corresponding induced representation for the Pauli amplitudes where the 
symbol V' is used to denote the usual nonrelativistic representation of the parity transformation, i.e. 

[P' ^^m] (PcP^) ^± (-P5,-p,). (242) 

Hence preserves parity, whereas T reverses parity, and, consequently, for a Salpeter amplitude with 

parity tt, i.e. V'J'm — positive energy Pauli amplitude has the same parity tt, whereas the negative 

energy amplitude has the opposite parity — tt: 

PVm=±^^m- (243) 

— The Salpeter amplitudes transform under rotations e 5*0(3), with rotation vector a; G H^, as 

K] iPi.Pv) = Su(S>Su(8>Su <Pii{RZ'Pi,RZ'Pn)- (244) 
In the Weyl representation we have 

Su= (j)^^ where u = exp (-i cr • u;) e SU{2), (245) 

with u (t{p) = <7{Ru) p) and we find the intertwining relations 

T^R^ = T±±± . (246) 

Here the induced representation 2?^^ of , which acts on the Pauli amplitudes (f^j , is exactly the usual nonrela- 
tivistic representation of the rotation R^ for a system of three spin-i fermions: 

[T^'r^ Vm] (p«> P'f) = u®u®u (p^.j{R-^p^,R-^prj). (247) 

Thus, to get the irreducible subspaces {J, Mj = — J, . . . , J} ofVji^, i.e. the Salpeter amplitudes with definite 
total spin J and 3-component Mj, the Pauli amplitudes ip^,^ have simply to be the usual eigenstates of the total 
angular momentum operator J = L -I- S as in a nonrelativistic system of three spin-i fermions. 

— As mentioned already, the three-quark embedding operators r='=±± explicitly break the flavor-S'C/(3) symmetry by 
the different quark masses m„ < and only a SU{2) (g) U{1) invariance remains. The operators of isospin and 
strangeness hence commute with the embedding maps T***, and accordingly the Pauli amplitudes ip^j are their 
eigenstates with quantum numbers T, Mt and S*. 

— The three-quark embedding operations T^^ ^ apparently are completely symmetric under arbitrary permutations 
a G S3 of quarks by their construction ( 236 ) , i.e. 

V„ r±±± = r±±± (248) 

where T)^ and T)'^ are the representations of the permutation o" G 6*3 on the Salpeter and Pauli amplitudes, 
respectively. This is a crucial point that permits to reduce the symmetry considerations of the Salpeter amplitudes 
to the embedded Pauli spinors. As the baryon Salpeter amplitude <P^.j must be totally antisymmetric, i.e. 

V, ^ij = sign{a) yaeSs, (249) 

the Pauli spinors themselves must have this symmetry: 

- sign((7) ^± Vae^s. (250) 
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In summary, the relativistic baryon Salpeter amplitude <Pfj m jTMtS* '^i^h specific quantum numbers J, tt, T, 
Mt and S* can be formed by embedding ordinary totally antisymmetric nonrelativistic baryon wave functions 

<PM J^MjTMtS'- 

+ ^ {Pi'Pv) V^M J-^MjTMtS'^P^'Pv ) (251) 

To define a basis for the totally antisymmetric three-quark Salpeter amplitude we can thus proceed in the same manner 
as in the nonrelativistic quark model, where the baryon wave functions fM m jTMtS' have the generic form 




j'MjTM^S'iPi^Pn) = 1 [V'Z(p«,P^)]k^ ® [Xsks ® ^ \ ®Ca> (252) 



1 






} (EC a] 









with 



— [V'l(P5' P'))]??,^ ^he momentum space wave function with total orbital angular momentum L, parity tt and permu- 
tational symmetry TZl G {S,A4s,Ma,-^}: 

— [xs]-fi the spin function of three Pauli spinors coupled to total spin 5* with permutational symmetry TZs G 

{S,Ms,Ma}; 

— [(f>Jj^ ]^ the flavor function with total isospin T, Ta-component Mt and strangeness S* which is of permutational 

symmetry TZp S {S,Ms,Ma,^}; 

— Ca the totally antisymmetric color-singlet state given by the Levi-Civita tensor: Ca = ^ciC2Ca |ci) ^ |c2) (8) jca). 

The momentum space wave function ipj^ and the spin function xs are coupled as usually to states of total angular 
momentum J, Mj according to 



{rLiPS,Pv)]TZ,^{xs]ns\^,, = E {LMl,SMs\JMj) rLM^Pi^Pv) XSMs (253) 

Ml, Ms 



with Clebsch-Gordan coefficients {LMl, SMs\JMj); the sum over the symmetries TZl, T^s and TZf m ( |252| ) is such 
that the combined momentum-, spin-, flavor wave function is totally symmetric, TZl Ei TZs 'S' TZp = S, and, finally, the 
baryon amplitude becomes totally antisymmetric with the totally antisymmetric color-singlet state Ca- 

Finally, let us discuss the implications of the present covariant approach with respect to the baryonic spectrum and 
the nonrelativistic quark model. According to the preceding discussion, the structure of the Salpeter amplitudes seems 
to be very similar to that usually considered in the nonrelativistic quark model: The Pauli amplitude embedded by the 
positive energy embedding map is of exactly the same structure as the usual nonrelativistic wave function for a 

given set of quantum numbers. But note the additional negative energy contribution to the Salpeter amplitude: For a 
specific parity tt of the baryon the embedding operator T brings also the nonrelativistic wave functions with the 
opposite ('wrong') parity — tt into play (due to the different behavior of T"*""*""*" and T under parity transformations). 
Thus, at a first glance, our approach seems to posses a larger number of states than the nonrelativistic approach. But 
recall that the negative mass (— M < 0) solutions of the Salpeter equation can be interpreted as the antibaryon states 
to baryons of just the opposite parity — tt, due to the CT'T-symmetry of the Salpeter equation. Exactly these negative 
mass states, which after CVT transformation become the baryon states (M > 0) with opposite parity, correspond to 
the additional states with the 'wrong' parity —tt. This feature becomes even more apparent, if we analyze the effect 
of the CT'T-transformation on the embedded Pauli wave functions. For the embedding maps we find the relations 

7V ® 7V ® 7V T***(P?, P.) = T TT^^(p^, p,) (254) 

such that the CVT transformation essentially switches the embedding maps of positive and negative energy and hence 
also the parity. Considering the decompositions of the negative mass solution ^^^^^ j,^^ g, 

^-M J^MjTMtS' " 7"+++ ip^M j^MjTMtS' +^ ^-M J-^MjTMtS' 

and of its related, CT^T-transformed positive mass solution of opposite parity (p^j j-^ m jTMtS* 

J-^MjTMtS* ^ ipl,j j~^MjTMtS* +^ ^M J^MjTMtS' 
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the CPT-transformation together with the property (254) of the embedding maps then yields the following relations 
for the corresponding Pauli amplitudes, 

'^tl J-^MjTMtS* = ~ '^-M J-^MjTMtS" 
'^M J^MjTMtS' = + "^-M J^MjTMtS" (^^'''^ 

which means that the CPT-transformation just interchanges the roles of both Pauli amplitudes. Thus we find in our 
present covariant approach exactly the same number of states as in the nonrelativistic quark model, a feature that in 
general can not be taken for granted in a relativistic approach: Consider e.g. the naive flavor-5J7(3) quark model. To 
explain the lowest lying multiplet it is assumed that the ground state orbital wave function is a totally symmetric 5- 
wave, the color state is completely antisymmetric and hence the spin-flavor state has to be totally symmetric, which in 
the nonrelativistic approximation restricts the possible multiplets to a flavor octet with spin- 1/2 and a flavor decuplet 
with spin-3/2. In a relativistic quark model, however, the number of spin-degrees of freedom is doubled for each quark, 
due to the presence of the lower components, which means that in the relativistic flavor-5'?7(3) model the number 
of possible symmetric spin-flavor multiplets is much higher than in the nonrelativistic approach, see e.g. |]8|, p^ , |3^ , 
Issll , in contrast to the experimental findings, which can be explained qualitatively by the naive nonrelativistic model 
very well. Our approach does not reveal this problem, owing to the Salpeter projector A and the CT'T-symmetry of 
the Salpeter equation, which circumvents such a proliferation of the number of states. Note that this feature of our 
model is a direct consequence of the instantaneous approximation or, more precisely, of the instantaneous ansatz for 
the genuine three-body (confinement) kernel. In view of the success of nonrelativistic quark models to account for the 
correct number of baryon excitations we in fact consider this to be one of the main empirical arguments to use this 
instantaneous ansatz. 



6 Summary and conclusion 

In this paper we presented how a relativistically covariant constituent quark model for baryons can be constructed 
within the general framework of quantum field theory. We started with the basic field theoretical quantities describing 
bound states of three fermions - the Bethe-Salpeter amplitudes and their adjoints - which form the residua at the 
bound-state poles of the six-point Green's function. The Bethe-Salpeter amplitudes, which might be considered as the 
covariant analogues of 'wave functions' in the ordinary nonrelativistic approach, obey a homogeneous eight-dimensional 
integral equation in momentum space - the so-called Bethe-Salpeter equation. In principle, this is the basic equation 
for the covariant description of bound states of three quarks in the framework of QCD, i.e. solving this equation for 
given single quark propagators and irreducible interaction kernels the discrete spectrum of baryons is then determined 
by the normalization condition. 

However, neither the full quark propagators nor the interaction vertices are reliably known functions in case of 
QCD such that reasonable phenomenological approximations for these basic ingredients of the Bethe-Salpeter equation 
were necessary. In order to remain as close as possible in contact with the features of the non-relativistic quark model 
we adopted the concept of constituent quark masses using free quark propagators with effective quark masses and 
the concept to describe the quark interactions by instantaneous, unretarded potentials. Although both replacements 
are chosen purely phenomenologically they are justified reasonably well by the apparent success of nonrelativistic 
potential models. As in the corresponding framework for mesons both assumptions then allowed a reduction of the 
full (eight-dimensional) Bethe-Salpeter equation to a reduced six-dimensional equation (Salpeter equation) in the 
case of instantaneous three-quark forces. In this case we obtained an equation for the reduced amplitudes (Salpeter 
amplitudes) with a structure quite similar to the ordinary Schrodinger equation and the normalization condition 
for the Bethe-Salpeter amplitudes reduced to the ordinary normalization condition inducing a positive definite 
scalar product. Complications arose, when two-particle interactions appeared, since these unconnected forces within 
the three-body system prevented a straightforward reduction as in the case of a pure three-body interaction alone. 
However, a reasonable treatment of these forces within the Salpeter framework is important since in quark models 
the three-body confinement forces are naturally supplemented by two-body residual interactions like the one-gluon- 
exchange or instanton-induced forces. We presented a method how in connection with the genuine instantaneous 
three-body kernel a reduction to a Salpeter equation of the same structure can nevertheless be achieved by deriving 
an effective instantaneous three-body kernel which parameterizes all effects of the two-body interactions. 

As a crucial property of the instantaneous approximation we found that it leads to a one-to-one correspondence 
with the states of the non-relativistic quark model, a fact which generally can not be taken for granted in relativistic 
approaches according to the doubling of the spin-degrees of freedom. In this respect the special projector structure 
of the Salpeter equation reduces the number of functions necessary to describe the bound state and thus circumvents 
a proliferation of the number of states: The Salpeter amplitudes, which still contain the full Dirac structure with 
positive and negative energy components, can be formed by an isometric embedding map of ordinary non-relativistic 
three-quark Pauli wave functions. 
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We found the appearance of the negative energy components to be related to the particle-antiparticle symmetry 
due to the CVT invariance: The spectrum of the Salpeter equation contains antiparticle solutions corresponding to 
particles with charge conjugated quantum numbers. This is a new feature of our Salpeter model for baryons and quite 
in contrast to ordinary nonrelativistic or relativized quark models. Solving the Salpeter equation for fixed spin J yields 
at the same time both the positive- and negative-parity bound-state spectrum and in particular positive and negative 
parity states are coupled in this way and are not independent as in nonrelativistic approaches. 

The fully relativistic kinematics and the formal covariance of our approach overcomes the old difficulties of non- 
relativistic approaches which in fact should be completely inadequate for small constituent quark masses. We expect 
the three-quark Salpeter equation to provide a more reasonable framework for quark models of baryons that should 
be superior to other treatments such as the nonrelativistic potential model or its simple so-called "relativized" 

extension |^ . In particular it offers the possibility to investigate the effects of the full Dirac structure of residual forces 
like the one-gluon-exchange or instanton-induced interaction and moreover it allows for the first time a reliable test 
of possible assumptions concerning the Dirac structure of three-body confining forces. In two subsequent papers 
we will therefore investigate explicit quark models based on the purely theoretical results of this paper and present 
concrete calculations of the complete non-strange and strange baryon spectrum up to 3 GeV. 
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A Appendix: Determination of the effective kernel V^f 

In this appendix we derive a prescription to construct the effective quasi potential V?f which has been introduced in 



sect. 14. According to eq. (157), V^f is defined by 

{Qm)a - A{gM)A, 



where on the left {Gm}a is given by eq. (159) which defines V^^ , i.e. 



{Qm)a = (GoAf) ~ i (GoA/) {Gm)a 



(258) 



(259) 



and on the right Qm is the solution of the integral equation (|142|) with the integral kernel Kf;j := Vn + K 



A3) ^17(2) 



'A/ J 



Gm ~ - i Gqm Kfj Qm- (260) 

Now the goal is to solve eq. (|25^ ) for V^f . For power counting purposes purpose we multiply the kernel Kf;^ by a 
parameter AG [0, 1], 

such that the Neumann series of Qm becomes a power series in A, and thus 

oo 

A{Qm)A = (GoAf ) + ^ A'^ A{GoM [-i K^] Gom ■■■ [-i ^m] Gom)A. 



(261) 



(262) 



fc=i 



k times 

The effective kernel V^^ becomes a function of A which is expanded into a Taylor series according to 



v; 



off 
M 



k T/cffC^) 
M 



k=l 



Inserting this into eq. (25£) the Neumann series of {Qm)a yields a multiple power series in A: 



{Qm)a = (Goa/) - i X! (^om) Vm'^ ^ (Gt 



OA// 



(263) 



(264) 



k=l 



oo oo 



+ ...+kr- 



{Gqm) 



{Gom) 



{Gom)- 
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Collecting all terms of equal power in the third term (which is of the order A-^), the multiple power series can be 
transformed into an ordinary series 



Ml 



(G, 



OA// 



OMI 



k=l 



(265) 



E E E (Gom) 

^■-2 fci, fc2, . . . , fcr < k 

ki + k2 + ■ ■ ■ + kr — k 



-iv; 



M 



(G, 



OM/ 



'iv; 



M 



(Go 



AI/ 



k=2 



Finally, we insert the resulting series ( |262| ) and (265) into eq. (258) and we arrive at 

oo oo / \ 

^ A'= (GoM> V^M (GoAf > = i E ^ ( [-i Go a/ ... [-i Kfj] Gom ) A 



k=l 



k=l 



k til 



E E E (GoA,) 

^==2 r-=2 ki,k2,...,kr <k 

ki + k2 + . . . + kr = k 



-iV, 



,fr(fci 



Af 



(G, 



OM) 



Af 



(Go 



M/ 



(266) 



which now enables us to solve for V^p{X) order-by-order by comparing the expansion coefficients of each power k of 
A. Thus, we find the following reduced terms which are irreducible with respect to (Go a/): 



— In lowest order, i.e. h — 1, we obtain the Born term V^f 



(1) 



(GoAf) Vm'^ "* (Gom) — ^{Gqm -^''m Gqm) A 



fr(k) 

and in fcth order, k > 2, we get V^j , determined by: 

(Gom) Vm'' (Gqm) = i ^ ( Go a/ [-i -S'm] Gqa/ 



-i Kf,] Gom ) A 



k til 



(267) 



(268) 



iE E 



(Go 



Ml 



'■-2 ki,k2, . . . ,kr < k 



Av; 



=ff(fei)' 



M 



(Go 



A// 



1 Vm 



(Go 



A// 



, ri2 , . ■ . , '>'r 
fcl -I- fc2 -I- . . . -I- fcr = fc 



Finally, we amputate the free Salpeter propagators (Go a/) using the Hamiltonian hoM with eq. (161) and thus, with 



the restriction (163) for V^f , we then can solve uniquely for V^f^'^K Consequently, we get the effective kernel as 



the following infinite sum of irreducible interaction terms V^^ 

oo 



M 
off(fe). 



where 



v; 



=ff(i) 



M 



^M 



hoM A ( GoAf GoM ) A hoM^ 



E < 



(k) 



k=l 



(269) 



(270) 



^M 



i hoM A ( GoAf (-i) 



G, 



OM 



(-i) 



T/(3) -I- 7?^^^ 



GoA/ ) A K 



OM 



k times 



E E 



A/ 



(G, 



OAf/ 



M 



(Go 



Af/ 



M 



L, . . . , . 

fcl + fc2 + . . . + fc,. = fc 



(271) 
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Finally, let us discuss how the instantaneous term i.e. that part of V^^^ that couples to the mixed energy 

components, enters in V^^'^^ As we mentioned already, this part of V^^^ appears solely in connection with More 

specifically 



1. In the Born term V^^^^ the isolated contribution of V^-' vanishes due to its instantaneity and its projector property 
AV^^^A = 0: 

(Gom V^'^ Gom) = (Com) vi'^ (Gom) = (Gqm) AV^'^A {Gom) = 0, (272) 

where we used (Gom) = (Gom) A = A{Gqm)\ 

2. For the same reason, reduced Feynman diagrams with more than two direct iterations of disappear: 

{Gom . . . Gom) V^^^ {Gom) V^'^ {Gom) V^^^ {Gom ■ • • Gom) = 0. (273) 

^ V ' 

=0 

Also, the ichIucx'cI irreducible kernel Vl(f does not contain two direct iterations of either, since such terms are 
reducible with respect to (Gom); because of 

^GoM • • • Gom Gom Gom ■ ■ ■ Gom 

= {Gom . ■ . Gom) V^^^ {Gom) V^^^ {Gom • • • Gom) (274) 

and thus are built by iterating two reduced Feynman diagrams of lower order. 

Therefore, we conclude that emerges in only such that is always directly attached to from the 
left and/or the right hand side. This means that in the Green's function {Qm)a (and even in Qm) at most two direct 

(3) (3) 

iterations of can occur. We want to remark here that this limitation of the number of direct iterations of V-^ 
offers an alternative counting scheme for the determination of via power series expansion, namely in powers of 

K^^ instead of powers of K^j = V^'^ + K^^ . 
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